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Abstract 


latex \noindent\vspace*{ 0.75cm} 

Oxygen steelmaking converter is the single most important unit process in a steel plant 
and is sometimes called the throbbing heart of a steel plant. Introduction of sublance tech- 
nology in the early eighties to analyze carbon and temperature in an oxygen steelmaking 
converter during the blow, approximately 3-4 minutes before the end blow, proved to be 
an important step in the control of oxygen steelmaking processes. Since the introduc- 
tion of sublance, the accuracy of end point prediction (hit rate) at most of steel plants has 
gradually increased from approximately 60% to 90%. 

In the present work several new models are developed by using genetic algorithm 
and multiple linear regression to predict end point carbon, temperature, phosphorus, man- 
ganese and dissolved oxygen from sublance data for 3 plants collected from literature. 
Other operational parameters like slag basicity, ore added during first blow, ore added dur- 
ing second blow, lance life, lance height, dolomite added during first blow etc. are also 
considered in predicting the end point values. 

Using the models developed in the present work, it is possible to predict end point 
carbon with an accuracy of cr=0.004% for Plant 1, a =0.005% for Plant 2, a =0.004% for 
Plant 3. Temperature can be predicted with <t = 5 °C for Plant 1, a = 8 °C for Plant 2, a 
=9 °C for Plant 3. Phosphorus can be predicted with a = 0.007% for Plant 1, a =0.001% 
for Plant 3. Manganese can be predicted with u =0.0065% for Plant 1, cr = 0.012% for 
Plant 3. Dissolved oxygen can be predicted with a = 34 ppm for Plant 1, a =140 ppm for 
Plant 2, a = 11 for Plant 3. 

The variation of capacity mass transfer coefficient K' is studied by using the dynamic 
models developed . The dependence of K' on gas flow rate is K' oc which implies 
that variation is slow and, during the blow, K' gradually decreases with time due to fall 
in rate of decarburization. It is observed that the variation is very slow with time. It is 
furthur observed that capacity mass transfer coefficient decreases with converter capacity. 
For example, K' for 100 1 converter is greater than that for 200 1 and for 300 1 converters. 
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Chapter 1 
Introduction 


Oxygen steel making process accounts for nearly 65% of the steel produced in the world. 
A typical diagram of the physical state of an oxygen steel making converter when a single 
hole lance is used, is shown in Fig 1.1. Most of the plants in world now use multi hole lance 
instead of a single hole lance.The pure oxygen jet emerging from the nozzle impinges on 
metal surface. A stream of metal droplets is ejected due to deflection of the jet in the 
impingment zone. The droplets eventually fall back to metal after their passage through 
slag or gas slag foam. The amount and size range of the droplets depend on lance nozzle 
design and lance height whereas the residence time of droplets in slag depends upon the 
properties of slag. The chemical reactions occur primarily in and around the jet impact 
zone and in the droplets falling through slag and rising gas bibbles. 

In the early years of its development (1952-1962), oxygen steel making process itself 
was given several names like LD (Linz-Donavitz), BOF (Basic Oxygen Furnace), BOP 
(Basic Oxygen Steel making), etc. These names were modified in early seventies when 
bottom blowing process was adopted in which pure oxygen gas was blown through tuyeres 
located at the bottom of the converter. The tuyeres needed cooling because of the high 
temperature generated at the tuyere tip. Depending upon the nature of coolant used (viz. 
fuel oil or natural gas), the new breed of processes were named as Q-BOP, LWS, etc. 
As a further modification of bottom blowing practice, combined top and bottom blowing 
processes were invented and today more than 30 commercial names exist for the variety 
of oxygen steel making process. For the sake of simplicity, however, they all are now 
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grouped under a single name of ‘Combined blown processes’. 

The specific combination of top and bottom blowing( i.e. the flow rate and nature of 
gas through top and/or bottom) adopted at a particular plant depends upon the starting 
composition of hot metal, capacity of converter, quality of lime, refractory life, desired 
end point composition of steel, economics of operation, etc. As a result, the steel making 
process at each steel plant has its own unique features. The blowing time may vary from 
15-22 minutes depending upon capacity, oxygen flow rate and nozzle design etc. The con- 
trol models developed at one plant are rarely applicable to the situation at another steel 
plant. In a more recent development (since mid eighties), a separate lance called “Sub- 
lance” is introduced into the converter to measure temperature and take a sample of metal, 
only few minutes before the end of the blow. Corrective actions are then taken on the 
basis of analysis of sample obtained by sublance measurement and this has significantly 
improved the accuracy of control of obtaining desired end point temperature and carbon 
of liquid steel. The fine tuning models based on sublance measurement, as also developed 
in this work, are called “sublance models”. 

The present work is organized into 7 chapters. Chapter 1 is the present one. Chapter 
2 briefly reviews the models published in literature as summarized else where [1], for 
end point carbon, temperature, dissolved oxygen, phosphorus and manganese, followed 
by the sublance models developed in the present work. Experimental data used in the 
present work are given in Chapter 3 and Chapter 4 discusses models for optimal control 
of end point windows. Chapter 5 presents results and discussion on the sublance models 
developed in this work and Chapter 6 presents results and discussion of the optimization 
models. Conclusions and suggestions for further work are presented in chapter 7. 
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Chapter 2 

Review of Models for End Point 
Carbon, and Temperature Prediction 
using Sublance Data 

2.1 Introduction 

Sublance is used in oxygen steelmaking to take a sample and also determine instanta- 
neously the steel composition 3-4 minutes before the end of blow. This analysis is used 
to determine the amount of oxygen to be blown in remaining part of the blow so as to 
arrive at a predefined carbon and temperature at the end. Thus, we are concerned with 
the modelling of trajectory of process during the last few minutes of converter operation. 
Metal sample is also taken at the end of blow to verify carbon, phosphorus, manganese, 
dissolved oxygen and temperature achieved. While the blow is in progress, adjustment of 
lance height can be done and judicious additions of iron ore, raw dolomite and solidified 
slag (return slag from converter) can be made to control the trajectory of carbon as well 
as temperature. Since the concentrations of carbon and oxygen are related to each other 
carbon may also predicted using the dissolved oxygen content of the bath and vice versa. 
Hence, the models for prediction of end point carbon, temperature, phosphorus, man- 
ganese and oxygen from sublance data play an important role in determining the amount 
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of oxygen to be blown during second blow period and thereby achieve the aim composition 
of steel. 

This chapter briefly reviews the models published in literature. Models developed in 
the present work are discussed in chapter 3. 


2.2 Review of models of decarburization, temperature and 

dissolved oxygen during end blow period 

One of the key objectives during end blow period (i.e. the period after sublance measure- 
ment) is the control of decarburization rate and a great deal of research has been done in 
this area. Several models, both static and dynamic, have been proposed which are suited 
to specific steel plant conditions. The fundamental basis of various models, adapted from 
an earlier work [1] are discussed in section 2.2 to section 2.4. 

2.2.1 Model by Lee and Rao 

Lee and Rao [2] have developed a mixed control model for decarburization in levitated 
droplets assuming gas phase mass transfer and dissociative adsorption of CO 2 (or interfa- 
cial chemical reaction) by studying the kinetics of decarburization in high carbon regime 
(0.91 1 - 2.479 wt % C) at 1973 K. The rate equation assuming mass transfer control in the 
gas, phase as 


dc 

dt 


1200 • A 
pV 


kg 

RTf 


In (1 + 


( 2 . 1 ) 


where Cc is wt % of carbon in the melt, A is interfacial area (cm^), p is density of the 
melt (^), V is Volume of the melt (cm®), kg is average mass transfer coefficient in gas 
film (^), Tf is average gas film temperature ( °K), is bulk pressure of C02. (atm). 
If interfacial chemical reaction control is assumed as a rate controlling step,then 


dCc 1200-^ (\ + KdfCs\ 

dt pV ^ ^[l + KfCs ) 


( 2 . 2 ) 
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where kR is rate constant (niol/m2/sec/atm), K is adsorption coefficient of Sulphur, 6^, is 
the fraction of the surface sites not available for the chemisorption of sulfur, / = ^ is 
activity coefficient of sulphur in the surface layer and in bulk metal, Cg is bulk 
concentration of sulphur (wt %), is interfacial pressure of CO 2 . 

On assuming mixed control of mass transfer and interfacial chemical reaction, the 
decarburization rate is 


dCc 

dt 


1200- A f kOkR(l + K0fCg) ^ 

[(k° + kR) + (k° + kRO^)KfCgJ ^ 


(2.3) 


where k° = 

The experimental results of Lee and Rao showed that 

•Liquid phase mass transfer is not a rate controlling step in the high carbon regime. 
The initial carbon concentration has no discernible effect on the decarburization kinetics. 

•Sulphur has a clear and reproducible retarding effect on decarburization. This effect 
is most pronounced in the range of 0 to 0.05 mass % sulphur. 

•The Pco 2 the gas mixture has a marked effect on the decarburization kinetics. 

•The flow rate of the gas mixture has a small but finite effect on the rate of decarbur- 
ization. 

In low carbon range (C< 0.4 mass %) equation assuming liquid phase mass transfer as 
control step is 


Pc = {n^c. - n*c) (2.4) 

where Rc is rate of carbon transfer (mol / sec), kc is mass transfer coefficient of carbon 
(m/sec), ric is concentration of carbon (b=bulk, i=interface ). 

Distin, Hallett and Richardson [3] have assumed in their study on the decarburization 
of liquid Fe-C alloys with O 2 and CO - CO 2 mixtures that the concentration of carbon 
at the surface becomes virtually zero at the instant of oxide formation. When the oxide 
appears on the surface Eq. (2.4) is reduced to 

kc.A ( h \ 
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(2.5) 


d[%C] _ M (oyr»^ 

dt V V J 

where A is surface area of droplet (m^), V is Volume of droplet (m^). 
For spherical droplets [2] 




( 2 . 6 ) 


where D is diffusivity of carbon, a is radius of the droplet. 

Lee and Rao calculated the kc for their experimental conditions as follows: 
T= 1973K 
( 7 * = 0 . 1 % 
a = 0.2875 cm 


D = 3.24 X 10-3 ^ 

sec 


— 


_ 5D _ 5 * 3 . 24 * 10 -^ _ 


a - 0.2875 = 56.348 * 10-^^ 


They found that mass transfer in the liquid phase is enhanced by about ten times due to 
electromagnetic levitation as compared to that at stationary condition. This enhancement 
is brought about as a consequence of motion in the liquid phase. A situation, similar to 
those for levitated droplets, may be assumed to prevail during decarburization in converter 
where droplets are ejected from jet impact zone at high velocity. 


2.2.2 Model by Chou et al. 

Chou, Pal and Reddy [4] have developed a general model for the decarburization process. 
If it is assumed that hot metal at 1400°C with a high carbon content is initially charged 
(for instance, liquid iron containing 4% carbon) in the converter with some scrap, then for 
the sake of simplicity bath temperature may be approximated to vary linearly with carbon 
content as follows: 


T = 1873 - 50[%C] 


(2.7) 



where T is the absolute temperature of the bath, [%C] is carbon concentration in the both 
expressed in wt. pet. 

In an actual converter, bath temperature can be predicted at each time step by dynamic 
heat and mass balance. When oxygen blows into the liquid metal, the carbon in hot metal 
will be oxidized by oxygen to form carbon dioxide. At an elevated temperature, the CO 2 
will react again with bath carbon according to the Boudouvard reaction 

CO 2 + [C] = 2CO (2.8) 

where [C] is the carbon dissolved in the hot metal. In accordance with the mass action 
law, we have 


K = 


p2 
^ CO 


(2.9) 


CO2 


where Pco, Pco 2 th® partial pressures of CO andC'02 respectively, Cc is activity of 
carbon in hot metal. K, the Equilibrium constant, is a function of temperature and can be 
derived in terms of standard Gibbs free energy data. 


P: = ea:p(15.30 - 16759/r) (2.10) 

Substituting Eq (2.7) into Eq (2.10) and combining with Eq (2.9) yield 


exp 


15.30 


16759 

1873 - 50 [%C] 


p2 

CO 

Pc02 ■ 


( 2 . 11 ) 


If the oxygen blowing into the converter bath is not pure but contains some inert 
gases such as nitrogen, argon etc. and the total pressure is assumed to be P, we would 
have 


Pco — P ' ^co 


( 2 . 12 ) 
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(2.13) 


Pco = P - X 


C02 


Pco = P-Xinr (2.14) 

Where Xco, Xco 2 , x/nr represent the mole fractions of CO, CO 2 and inert gas in the 
decarburization zone, respectively, and 

x^2 Xjjij' ~ 1 (2.15) 

On the other hand, the activity of carbon can be expressed in terms of the products of 
[%C] and activity coefficient, fc i.e. 


ac = [%C]fc (2.16) 

Substituting Eqs (2.12) to (2.16) into Eq (2.11) generate as 


exp 


15.30- 


16759 


P-xl 


(2.17) 


1873 - 50 [%C] J (1 - Xinr - Xoo) [%C] ■ /c 
Let two other quantities be defined as Xco and Xco 2 (which are relative contents of CO 
and CO2 in CO + CO2) 

Xco=Xco/(Xco+Xco2 )~X(^/( 1 ~X/jjf.) 

Xco2~^co2 ! (^co"*"^co2 )— ^C02(( 1 "^/nr ) 

Substituting XcoandXco 2 into Eq (2.17) we obtain 


exp 


15.30 


16759 

1873- 50 [%C] 


P • Xj (1 - X,..) 

(l-Xco)[%C]-/c 


(2.18) 


The above relation gives the dependence of relative content of the carbon monoxide 
on the carbon concentration in the hot metal during the BOP operation. Based on this 
solution the relative mole fraction of CO can be calculated if carbon content and fc are 
known. 
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Bensel et al [5] have used the following expression for f^. 

fc = a[%C] (2.19) 

where a = 0.167 when [% C] > 1%, a = 0.2 when [% C] < 1% 

These are approximation formulae and they completely neglect the effect of temper- 
ature on fc- Rist and Chipman [6] have given a Equation, which explains the effect of 
temperature on activity coefficient of carbon 

\ogfc = 0.1666 [% C] - 0.01585 [% C]^ + 9.9613 X lO'^ [% C]3.(T-273) 
+3.0246X10-^[% C]^-(T- 273) 

The decarburization rate depends on the oxygen blow rate into the furnace. If the rate 
of oxygen blow is (^) the decarburization rate then from mass balance 
considerations 

(1 ~ ^Inr) {^C02 + 2^co) 

Where the factor (1 - X/nr) reflects the purity of oxygen. The negative sign stems 
from the decreasing carbon content during decarburization. It is implicitly assumed that 
oxidation of other elements like P, Mn, Fe, S is negligible in high carbon regime and this 
is approximately so. 

Expressing in terms of [%C] and the flow rate Qq (Nm^ /ton) the above equation 
becomes 


d[%C] ^ -Qo (1 - Xinr) 
dt 1.87 ■ 10-5 (1 _ o.5Xco) Wt ^ ’ 

where Wj is bath weight in gms, [%C] is weight percent of carbon. 

When carbon content is less then 0.3 , the decarburization will be controlled by mass 
transfer of carbon in the melt instead of oxygen and it can be expressed as 

^ = ( 2 . 21 ) 

where, C**, (7* are bulk and interfacial concentration of carbon, mol/cc, kd is mass 
transfer coefficient (cm/sec), A is bath / oxygen interface area (cm^). 

Since the reaction between carbon and oxygen at high temperature is very fast, C® is 
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negligible (except towards end of the blow) and hence Eq (2.21) can be written as: 




or 


where V= melt volume 


d[%C] 

dt 


kfi • A 
V 


[%C] 


( 2 . 22 ) 


2.2.3 Exponential Model 

2.2.3.1 Model by Bessho et al. 

Bessho et al.[7] have developed a dynamic model, based on exponential dependence of 
decarburization rate on carbon content . This model has been applied to estimate carbon 
content and temperature of Q -BOP heats after sublance measurement. 

Exponential model is expressed as 


where 


dC 

d02 


= ^-t-rexp {(l>C) 


(2.23) 


r = -^ exp 




Then Eq (2.28) becomes 


dC 

d02 




1 — exp 


' Co -C 
Cp-Co 


(2.24) 


where ^ is theoretical maximum decarburization rate, Co=Critical carbon content below 
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which decarburization does not substantially occur and is found to be 0.01% for Q- BOP, 
Cp = a parameter derived statistically based on operational variables. 

Value of can be calculated for a converter with oxygen flow rate Qo = 600 Nm3/min 

and weight of bath Wj = 300 tons as follows: 

600 1 1 9.9 

f = iSSS = -IWr^Imm 

By integrating the Equation (2.24) 



where Csz. is carbon content determined by sensor lance, 

C/a is the aimed carbon content, 

A (92 is the amount of oxygen required (per ton of charge) to reduce the carbon content 
from C 5 Z, to C/a . 

End point carbon content at any arbitrary time after the measurement by the sensor 
lance can be calculated, by simplifying the above equation to 


C/a = (C'a-C'p) In 




(2.26) 


where. 


^A02 


~ CsL 
Co 


where AC 2 (Nm^ /ton of charge ) is the amount of oxygen blown after the measurement 
by sensor lance. The value of Cp in the above equations is estimated from multiple linear 
regression (MLR) analysis which involves the carbon content and temperature of the bath 
monitored by the sensor lance.The MLR equation to estimate Cp can be expressed as 


m 

Cp — (Xq ^ * x^i 

1=1 


(2.27) 
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where aoA = regression constant and regression coefficients, Xj = operational variables as 
shown below. 

The operational variables used by Bessho et al. For regression analysis are 
Xi meaning 

1 real carbon content at blowing end 

2 real bath temperature at blowing end 

3 hot ratio (=wt of bath /wt. of scrap ) 

4 carbon content by sensor lance 

5 bath temperature by sensor lance 

6 amount of lime used after one point check by sensor lance. 

7 total amount of lime and dolomitic lime used 

8 calculated basicity 

Bessho et al. used a semi empirical equation to estimate the temperature rise. 


Tf — Tsl + P • ^0-2 4 - Q 


1 


CsL 


+ R 


(2.2S) 


where T/ , C/ is the temperature and carbon content of the bath at blowing end 
respectively, 

Tsz, ,CsL is the temperature and carbon content of the bath measured by the sensor lance, 
AO 2 is amount of blowing after sensing ( Nm^/ton of charge), 

P, Q, R = Parameters calculated from multiple regression analysis by using actual 
operational data such as lance height, slag 'volume , etc. 


2.2.3.2 Model by Byun et al. 

Byun et al. [8] employed a simplified form of exponential model Eq. (2.24), as: 

- = a [1 - exp -/? (c - Co)] (2.29) 

Wst = Weight of molten steel (kg). Here the term/5 is same as l/(Cp -Co) and a is same as 
^ in Eq (2.24).It is implied that 13 represents combined effect of Cp and Co on 
decarburization. 


12 



Integrating Eq (2.29) from sublance carbon content Cj to aimed carbon Cg gives the 
oxygen consumption from the in - blow sublance measurement to the blow end 


^02, cal — ~ [ ^^2 = — [ 

JCs JCi 

Eq (2.30) can be simplified to 


W. 


St 


cs [1 - exp -P (c - Co)] 


d[C] 


(2.30) 


AO2, caJ — 


Wji 

OtP 


/3(C'e-C,)-ln 


l-exp-/?(Ce-C,) 
1 - exp -P {Cs - Ce) 


(2.31) 


The effect of other operational variables on oxygen consumption can be incorporated 
through multiple linear regression as 


%(X,-X,)+7o 

(2.32) 

where = regression coefficient, 70 = regression Constant, X, = operational variables, 
Xj= Standard operational variables. 

2.2.3.3 Temperature calculation model 

Temperature at the end of blow can be calculated by considering the heats evolved in 
various reactions 


^02, cal = 


w, 


st 


aP 


P {Ce - Cs) - In 


1 - exp -P {Ce - Cs) 
[1 - exp -P {Cs - Ce)J 


+ E 


[C] + \02(g) = CO(g) 

AHi 

(2.33) 

CO(g) + i02(p) = CO 2 

AH 2 

(2.34) 

fe+iOj(s) = (FeO) 

Ai?3 

(2.35) 
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(2.36) 


^02 = [O] AH, 

Vhere, AHi is the heat produced from each reaction (Kcal / Nm® - O2) 
Over all heat balance (excluding heat losses) can be written as 


CpWsfdT = AH\dOc + AH2dOco + AH^dOpe + AH,dOdts (2.37) 

^here Cp is specific heat capacity of molten steel (Kcal / kg °C),Oc = oxygen combined 
dth carbon to form CO (Nm® - O2), Oco = oxygen combined with CO to form C02 (Nm^ 
O2), 0^6 = oxygen combined with Fe to form FeO (Nm^-02),0dts = oxygen dissolved 
i liquid steel (Nm^-02 ). 

In Eq(2.35), the heat of Fe oxidation and the heat produced by oxidation of CO to CO2 
e not greatly different at steel making temperature, i.e. AH 2 « AHz. In the later stage 
' EOF blowing process , it is generally accepted that the product of dissolved carbon and 
tygen in liquid steel is constant. 

[C]x[0] = K; 

Therefore, the amount of dissolved oxygen can be expressed as a function of carbon 
ntent in steel. 


dOdis • 


22.4 • Wst 
3200 


22.4. K', 
3200 



(2.38) 


Assuming that the carbon removal from steel melt mainly occurs by Eq(2.33), the 
„ 'gen consumption for decarburization can be obtained as 


dO, = - [C] = -K,d[C\ (2.39) 

f an assumption is made that all lance oxygen is consumed by reactions (2.33) to 

0 , 

d02 = dOc + dOco + dOrfts 
5 term dO/re + dOco is expressed as follows. 
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dOpe + dOco — d02 — {dOc + dOd^s) — d02 + i^ic?[c] + i^4d[c]/[c]^ 


(2.40) 


If iJ j are assumed to be independent of temperature , the following differential equa- 
tion can be obtained from Eq(2.37) through (2.40). 

dT = ( (di + (da / [Cf) d[C] + dadOa 

where di= ( AH 2 - AFi)Ki/(Cp W,* ), da = ( AH 2 - AF 4 ) K 4 /(Cp W,t),d 3 = AH 2 / (Cp 
. Wji). Integration of the above equation from in blow measurement(rs) to the end blow 
(Te) gives a linear function of temperature rise. 

AT,ai = di (C, - Ce) + <^2 (^ - ^) + daAOa (2.41) 

where O 2 is total oxygen blow between these measurements. 

In order to fit this model better to the actual temperature data, parameters dj in Eq(2.41) 
are modified by the multiple regression analysis as discussed in the decarburization mod- 
el. In the present work operational variables other than the amount of oxygen and carbon 
content were also included in this model such as the steel temperature at sublance mea- 
surement, slag volume, HMR, lining life, lance height and the amount of iron ore used. 
Therefore, the equation for the temperature increase can be rewritten as; 

ATactual = €1 [Cs — Ce) + €2^ _ q + €3A02 + Xt {^t ~ ^i) + XO (2.42) 

where Xt, is regression coefficients, xo is constant, Zt is Operational variables and^ is 
standard operational variables (mean values ) 

2.2.3.4 Feedback model 

There will always be deviation of oxygen quantity and temperature rise for a heat which 
are not accounted for by the explanatory relationships such as Eq. (2.32), (2.41) and 
thus some part of those deviations will remain unpredictable. In order to minimize the 
error various steps can be taken. The model can be modified to include the perturbing 
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factor which may be explained by the selected model parameter. The model parameter 
for upcoming heat can be adjusted by a time series analysis of the posterior calculation of 
actual turn down data of preceding heats. In other words the new model parameters for 
upcoming heat can be expressed as a function of old parameter of preceding heats. 

One of the most important concerns in developing feedback model is how to select the 
adapted model parameters. In this work, various model parameters, 0, for decarburization 
model and 4>i for temperature calculation model were employed and the general form of 
equations used are 


AO2 — di 



'1 - exp 03 (Cs - Co)' 

.0203 ^ 

1 - exp 03 {Ce - Co) 


t 


(2.43) 


AT = <l>i 


Cl (C's ~ C'e) + €2 + <j>2^02 + ^Xt (-^i ~ 


(2.44) 


In addition to above, the following time - series forecasting methods can be tested 
to obtain appropriate feed back model. 


Method 1 ; Single moving average This technique consists of taking a set of observed 
values, finding the average of those values, then using that average as the forecast for the 
next period. 

Ft+i=(Xt+Xt_iH + X t-N-i ) / N 

where F t+i = the forecast value for the next period, Xt = the most recent observation, N 
= the number of historical data. 


Method 2 : Exponential smoothing method. This forecasting is based on weighting 
the most recent observation with a weight of value,K* whose value is between 0 and 1 and 
weighting the most recent forecast with a weight of value, (1 - K*). 

Ft+i=K*Xt+(l-K*)Ft 

where K* = Smoothing factor 
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Method 3 : ARMA method In this model, the series to be forecast is expressed as a 
function of both previous values of the series ( anti regressive term ) and previous error 
values from forecasting ( moving average term ) . ARMA (2, 1) model used in Byun et al. 
model was 

Ft+i='^iXt+A2Xt_i-/iie( 

where Ai, are coefficients, e* is the most recent error from forecasting . 

Of the above three methods, the exponential smoothing method has been reported to 
give the best result with a K* value of 0.1. The procedure for On-line application of the 
model is as follows: 

The in- blow sublance measurement is made when 80% of total oxygen amount, pre- 
sumed by the static, model has been blown. At this point, the decarburization model starts 
a calculation of the oxygen amount to be blown for the aimed carbon content. The temper- 
ature increment in this dynamic period is then calculated based on the amount of oxygen 
to be blown. If the calculated end- point temperature is higher than the aimed one, the 
amount of coolant to be added is determined. If the calculated temperature is lower than 
the aimed one, the heat has to be over blown. The bath temperature and the carbon content 
in steel at turn down are used for adapting the model parameters for the next heat. 

2.2.4 Model by Szegedi et al. 

Szegedi et al [9] have developed a model for the adjustment of carbon content of steel by 
dividing the decarburization process into 3 periods, namely 
I . Thermodynamic period 
n . Reaction kinetic period 
III . Chemical period 

as shown in Fig. 1.1. V^c in the figure represents decarburization rate given by 


The formulae to calculate the decarburization rate and time taken for the blow for each 
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[%C] 

Figure 2.1: Schematic representation of decarburization process. 


period are as follows. 

The decarburization rate for the period I is obtained on the basis of operational data as 
Vc./=0. 1071. i.[-0. 1421. CB-0.7366.SiB - 0.0425.i + 1.4987 ] 
where Vcj is rate of carbon oxidation(%C/min), i intensity of blowing (m3 o2/t s- 
teel.min), Cb carbon content of charge (%), Sia is silicon content of charge (%). 

Time taken for computing the period I is given by 

t/=100-SiB[i'(2.03TCB+7'SiB + 0.106 •! - 6)] min 
Similarly for the reaction kinetic period 11 

Vc,//=0.1071-z 

and time necessary for second period is 

t//=9.375(CB-0.046)/(-1.942-Si.A/(i-B)) 
where A= Cb +5.1728.SiB + 0.2981. i -10.52426 


and B = 0.4217 - 0.1426. Cb - 0.4918 - SIb -0.00744.i 

The third period is chemical period for which the rate of carbon oxidation is 

Vc,j//=12.89[%C]-exp(-6590.73/T) 

The regression analysis of operational data has shown that 

_ exp 2 l 4 ^zf^ 0.208i+0.046 

tjll - 1289 [%Cl > 

The activation energy of the reaction is 
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G^=-R-6590.73 J/mole. 

Total time for blowing t = t/ + 1// + 1/// 

By equating the decarburization rates periods II and III (i.e. Ycji = ^c,iii ) they have 
derived a formula to calculate the critical carbon constant [Cjt]. 

Cfc=0.208.(i)+0.046 

They have experimented this model in the range of 0.01 -4.0 C wt % .and have shown 
that the adjustment and the more exact control of carbon content ensure a significant saving 
of time and material costs. 


2.2.5 Models Based on Exhaust Gas Composition 

2.2.5.1 Model by Fukumi et al 

Fukumi et al [10] have developed a refining control model based on exhaust gas informa- 
tion measured by mass spectrometer. This model consists of calculation models that use 
sublance information to estimate various values as follows. 

•Calculation and estimation of the oxidation degree ( W02) 

•Estimation of (T.Fe) 

•Estimation of decarburization 
•Estimation of temperature in steel 
•Estimation of oxygen content in steel. 

2.2.5.2 Calculation and estimation of oxidation degree (WO 2 ) 

WO 2 is defined as Eq(2.45) based on oxygen balance calculated by refining conditions (top 
and bottom gas flow rate, amount of flux, etc.) and exhaust gas information (analysis, flow 
rate) 


VF02 = 


/ (fi(Q,F)- jo(V,A))c£^IWM 


(2.45) 


where fi is input oxygen, fo is output oxygen, WO 2 is oxidation degree, Q is top blow 
oxygen, F is amount of coolant, V is 02,C0, CO 2 volume in exhaust gas, A is sucked air 
or slopping gas volume. 
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WO 2 indicates oxygen in slag and steel, and this value is the most important parameter 
in this model. The volume of air sucked from the converter mouth and slopping gas volume 
are evaluated with reference to bottom blowing gas, balance of air, Ar and N 2 . Analyzed 
values of exhaust gas have a time lag behind the actual reactions because of gas transport 
time. At the final stage of blowing when W02 is rapidly increased, this time lag must 
not be neglected. Thus, reacted volume can be calculated at real time by the estimation of 
W02 from Eq(2.46) which is calculated from actual data. 

dW O 2 / d02 = a [1 - exp W O 2 /P] (2.46) 

where WO 2 is degree oxidation, O 2 is input oxygen, a constant, /3 is parameter. 

2.2.5.3 Estimation of total iron 

Since WO 2 indicates oxygen in slag and steel , WO 2 has a close relationship with total 
iron (T.Fe) in slag. In other words, the change of (T.Fe) can be known from the change of 
WO 2 . Therefore, from (T.Fe) at the first sublance and the increasing amount of WO 2 after 
the first sublance, the value of (T. Fe) at the final stage can be estimated in real time from 
the equations given below. 

(T.Fe)s=f(QB,F 02 ,Cs) 

Ws(t)=Q!W 02 (t)+PF,<,/{ l-^(T.Fe)s/100 } 
(T.Fe)(t)=(T.Fe)s+ 7 W 02 (t)AVs~t 

where, (T.Fe)s is (T. Fe) at in - blow sublance, QB is bottom blowing gas flow rate, Cs 
is [C] in steel at in - blow sublance, WO 2 is oxidation degree, Wgo is Slag volume except 
(FeO), (MnO) , (P 2 O 5 ), Ws is slag volume, t is time in sec, a, 7 are parameters. 

2.2.5.4 Estimation and prediction of the amount of decarburization 

The amount of decarburization utilized for the estimation of steel temperature and pre- 
diction of slopping is calculated based on carbon balance from refining conditions and 
exhaust gas information, from the equation 

AC=[/ /(ho(V.A)-hi(F)dt]AVM 
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where, ho is Output carbon, hi= Input carbon, AC is Amount of decarburization, V is 
CO, CO2 volume in exhaust gas, A is sucked air volume or slopping gas volume, F is 
amount of coolant and flux. 

At the final stage of blowing , the estimation of decarburization , amount that contains 
time lag is calculated according to the efficiency of oxygen for decarburization. 

2.2.5.5 Estimation of steel temperature 

The heating of the steel at the final stage of blowing has a close relationship with the 
generated amount of (T.Fe) namely the increased amount of WO2 and amount of decar- 
burization. 

Therefore, heating after sublance is estimated by 

dT/dOa =a . dW02 / 602-^ . dC / d02-7 WsWM + 5 
where T is temperature in steel, O2 is input oxygen (gas + coojlant), WO2 is oxidation 
degree, C is carbon content in steel, Ws is slag volume except (FeO) , (MnO) , (P2O5), 
WM is Weight of steel, a, /?, 7, 5 are parameters. 

2.2.5.6 Estimation of dissolved oxygen content in steel 

Oxygen content in steel is estimated from the oxygen balance Eq (2.47) calculated from 
(T.Fe) which is evaluated from exhaust gas information and the steel temperature. 


log[0] = log [a (T.Fe)] -f ft/T + 7 (2.47) 

where, [0]= Oxygen in steel, T is temperature in steel, a, P, 7 are parameters. 

2.2.6 Model by Iwamura et al 

Iwarmura et al.[l 1] have developed a new endpoint control system of EOF with an auto- 
tuning function for the model parameter based on the measured exhaust gas composition. 
This system can accurately estimate the steel temperature and the carbon content at the 
blow-end through one campaign without changing of the model parameter. The dynamic 
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control model consist of an oxygen balance equation and heat balance equation. The 
amount of oxygen required in a heat is expressed by the sum of the amount of oxygen 
based on the change of carbon. 


- = ao (C, - Ce) + ai In ^ + a2 + 13 k {X, - X^) + FB, (2.48) 

where ao,ai are coefficients of oxygen balance, Cs, Cg are carbon contents at sublance 
measurement and endpoint, FO2 is the amount of oxygen after sublance measurement, ri 
is coefficient of oxygen in fluxes, Wsub is amount of flux, O2 is the amount of oxygen 
fluxes, W ST is weight of molten steel, X^ is factor i of operating condition, Xj is mean 
value of Xi, k is coefficient of statistical model, FBc is adjustment of oxygen balance 
from the feedback. 

In the heat balance model, the fundamental equation has been changed from the carbon 
content to the oxygen basis . In situations where the carbon content of the metal at the 
endpoint has little variation, temperature variation has stronger relation with the amount 
of oxygen than the change of carbon content. The heat balance is represented by the 
following equation. 


Te - r, = bo + 7 E WsTsub + E (2.49) 

where Ts,Te. is temperature at sublance measurement and at the endpoint, bo, b2is co- 
efficient of heat balance, 7 coefficient of cooling capacity, k is coefficient of statistical 
model, Tsub^s cooling capacity of fluxes, FBt is adjustment of heat balance from the 
feedback. 

2.2.6.1 Auto-Parameter -Hining function 

Usually the decarburization efficiency is changed during blowing for each heat. This fact 
shows that in the oxygen balance equation with a fixed parameter it is difficult to obtain a 
highly accurate estimate of steel carbon content and temperature at endpoint through one 
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campaign . To solve this problem, Iwamura et al. have added an auto-parameter -tuning 
function of parameters in oxygen balance Eq.(2.48) using exhaust gas information from a 
mass spectrometer. The parameters of Oq and are adjusted according to the following 
Equations. 

First of all, Eq(2.50) is altered to the following Equation . 


—W ST 


dt 



(2.50) 


The decarburization efficiency is expressed with the parameter of oq and and the 
estimated steel carbon content . On the other hand, the actual decarburization efficiency is 
calculated by 


a = G 


WsT A C^ 
Ao^ 


(2.51) 


where G is gain, a is decarburization efficiency, i is sampling number from sublance 
measuring time to blow-end. 

ACi is the change of steel carbon content and AOj is the amount of oxygen during the 
sampling period At as shown below. . 

AQ=Q-Q_i 


Aot = 



(2.52) 


The amount of decarburization in the exhaust gas during (t is calculated from the 
following expression using the exhaust gas information. 


AWr 


Jf 


t'+At {CO (r) -b CO 2 (r)} Q (r) 12 
3600 • 22.4 • 1000 


dr 


From Eq. (2.51) and .(2.53) we can obtain 


(2.53) 
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(2.54) 


Ao, 


when C is infinity, the term of 1/C in the Eq.(2.50) can be neglected and the parameter 
of ao is expressed fi'om Eq.(2.50) and (2.51) can be shown as, 


ao = — (2.55) 

oci 

In on-line operation, the parameter of oq is decided by the average of the actual decar- 
burization efficiency during the maximum decarburization period before sublance measur- 
ing time as shown in the following equation 


CLq 


1 


(2.56) 


where, oimax= decarburization efficiency in maximum decarburization period, m is 
constant of sampling number. 

On real time after sublance measurement, the parameter of ai in Eq.(2.50) is calculated 


as. 


ai 


(l-ao-aQ ^ 

a* 


(2.57) 


The parameter of ai changes in proportion to the decrease of the decarburization effi- 
ciency and this function can accurately estimate the steel carbon content. 


2.2.6.2 Feedback calculation model 

The feed back calculation model has been developed in order to adjust the errors in the 
dynamic control model. The amount of feedback is calculated by considering the trend 
in the errors of the dynamic control model in past some heats measured by sublance at 
blowing and endpoint. The rule of feedback is classified with three cases by the trend of 
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the error as shown below. 


Case 

The trend of the error of the model 

The amount of feedback 

1 

Mdlerror(i- 1 )-Mdlerror(i)>0 

FB(i) = FB(i- 

1 )+G F Berror(i)+error(i- 1 )/2 

2 

Mdlerror(i-1)-Mdlerror(i)<0 and 
lerror(i)l-lerror(i- 1 )>H 

FB(i) = FB(i-l)+GpB-eiTor(i) 

3 

Mdlerror(i-1)-Mdlerror(i<0) and 
lerror(i)l-lerror(i- 1 )l <H 

FB(i) = FB(i-l) 


where, FB(i) is feedback value, Mdlerror(i) is error of the model, error(i) is error of the 
feedback model, is feedback gain, H is on-detective value. 

In the oxygen balance equation, the feedback value is adjusted according to the con- 
dition of the sublance measured carbon content at end point in addition to the above rule 
of the feedback. This is because the sensitivity of the carbon content for the amount of 
oxygen varies with the sublance measured carbon content in the blowing. 

The relationship between carbon content and amount of oxygen can be defined in the 
following equation 


Gbp = Si log (2.58) 

where si, $2 are coefficients. 

The sensitivity of the carbon content to the amount of oxygen is related to the following 
equation. 


d A 02 

dCsp 


1 

— exp 
5i 


Cep ~ S2 
S2 


(2.59) 


where e is the sensitivity 

The actual feedback value is adjusted with the sensitivity as shown in the following 
equation. 

FB’(i)=h-FB(i) 

where 

h = -^ 

^base 
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and 


^ — — exp 

cluep Si ^ S2 


2.2.7 Droplet Decarburization Model 

This model was developed by Divyesh Dixit [12] .The metal droplets are ejected from jet 
impact zone and eventually return to metal via slag phase. The dissolved carbon in these 
droplets reacts with FeO present in the slag( or FeO formed on the surface of droplets) 
according to the following reaction 

[C] + Feo = {CO)g + Fe 

If the droplets spend on average r time (residence time in slag foam) then 


In 


' C^-Ce ' 

Ci - Ce. 


KaA 

-nr-T 


(2.60) 


Where is the carbon concentration of droplet at the time of falling in metal, is 
droplet carbon concentration at the time of ejection of droplet, Ce is equilibrium carbon 
concentration, A is surface area of the droplet, V is volume of the droplet and Ka is the 
mass transfer coefficient within the droplet. 

From the above equation. 



1 — exp 



where is carbon concentration of droplet at time of falling in metal 
For the droplets joining the metal bath between time t and 1 4- At, 


(2.61) 


d = 


where the concentration of vVcarbon in bulk metal at time f — r. 

Hence Eq. (2.61) can be written as 


cf = a 


1 o 

1 - exp — — T 


I K^A 

+ Cl exp — —T 


(2.62) 


where 
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If iron conversion is defined ad Kg/s of the droplet thrown in to the slag, then the mass 
of droplets produced between time t and / + Af is Mg.dt. From carbon mass balance, 


WbC^ + Mg ■ dtCf = (Wfc + Mg • dt) C; 


(2.63) 


this can be simplified to 


\ C§-Ce ] _ Mff[l-exp-:^r] 

[ci-Ce\ Wb + Mg-exp-^r ^ 

where, Wb is bath weight, Mg.dt is droplet fall in metal,C^is bath concentration before 
fall, C^is bath concentration after fall. 


2.2.8 Quadratic Model 

This model, proposed by Besso et al [7], postulates that at a later stage of blowing oxygen 
gas is consumed in three different ways: (1) formation of CO, (2) formation of FeO and 
(3) formation of other oxides and stray oxygen. Assuming that consumption of oxygen 
for (3) is very small, we can write 

d02 d02p,o + cJ<^ 2 co (2-65) 

where d02 is the total oxygen, d02Feo amount of oxygen consumed for FeO forma- 

tion and d02co the amount of oxygen consumed for CO formation. 

The concentration of FeO in the slag is assumed to be inversely proportional to the 
concentration of carbon in the bath i.e., 

(Fe(7)= ^-t-7 

c 

where and 7 are constants and statistically determined. 
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Suppose percentage of oxygen which forms FeO is represented as T]Feo (oxygen effi- 
ciency for FeO formation), and the percentage of oxygen which forms CO is represented 
as T]co (oxygen efficiency for CO formation). Both rjpeo and rjco are related by a param- 
eter a as 


VFeO = Oi (100 — rjco) (2.66) 

where 0 < a < 1. The parameter a in eq. (2.66) is introduced to take care of the approx- 
imation in eq. (2.65). The value of this parameter can be calculated from experimantal 
data. Theoritical amount of oxygen required for CO and FeO formation can be calculated 
as given below. 


dO^co = -d [^c] * 1000000 * 22.4 * 10-^ * ^ 
= —9.33 * Wm * dc 

where Wm is the weight of metal in tons, and 


1000000 * 22 . 4 * 10 -^ 

72 


k^W.d[(FeO}] 


where Ws is the weight of slag in tons. Substituting (FeO) in the above equation, we get 

1.556 * p 


dO- 


2FeO 


-dc 


a * 


Substituting d02co (2.65) we get 

'1.556 


d02 = — 


* a 


-b 9.33W„ 


dc 


(2.67) 


since the rate of decarburization — ^ in the eq. (2.67) is expressed as quadratic funciton 
of carbon content in the bath, this model is denoted as quadratic model. Upon integration 
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in the limits csl ans cpa we get 


AO 2 = 9.33W„, (C5L - cpa) + 1.556W^,- (- (2.68) 

\CsL ^FaJ 

where AO 2 is total oxygen blown, csl is sublance carbon and cpa is end point carbon. 

One of the limitations of this model is that it neglects the oxygen required for oxidation 
of P, Mn etc, also neglects the oxygen for oxidation of CO to C'02. The percentage of 
CO2 in waste gas can be as high as 15-25% towards end of the blow. 
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Chapter 3 

Models Developed for End Point 
Carbon, Temperature, Phosphorus, 
Manganese and Dissolved Oxygen 


3.1 Introduction 

The flux of carbon in metal, assuming mass transport of carbon in metal to be the rate 
controlling step, is given by 


or 


— — —k—ic'’ — d) 

dt- V'" 




(3.1) 


where A is the area of interface, V is the volume of bulk metal, t is time in seconds, d 
is concentration of carbon in the bulk, c* is the concentration of carbon at the interface 
and K' is capacity mass transfer coefficient (or rate constant). Capacity mass transfer 
coefficient K' depends essentially on stirring energy imparted to metal by impinging jets 
and stirring energy due to gas evolution within bath and bottom stirring. In a small time 
interval (for constant lance height and gas evolution rate) if k is assumed to be constant 
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then eq. (3.1) can be integrated to give 

In = -K'-^ (3.2) 

where cj is end point carbon, Co is sublance carbon and Cg is the equilibrium concentration. 

It may also be noted that, besides gas flow rate, K' may also dependent upon opera- 
tional parameters such as basicity, hot ratio, ore added in first blow (orel), dolomite added 
during first blow (dolol), slag added during first blow (slagl), lance height during second 
blow (hlans 2 ), ore added during second blow (ore 2 ), dolomite added during second blow 
(dolo2) and slag added during second blow (slag2) etc. For small concentrations of c* in 
comparison with eq. (3.1) can be simplified to 


= -K'c 


It can be shown that stirring energy density (e, is directly proportional to gas flow 
rate[13] through metal and the capacity mass transfer coefficient is proportional to 


€ a (gas flow rate) 

Thus Ki can be written as 

r dc V 

X'oc 

where — ^^Siis rate of decarburization i.e., gas (CO) evolution rate and Qe/f is effective 
volume flow rate of gas from top and bottom (see appendix 1 ) and rj is exponential 
factor. The magnitude of exponent 77 can either be assumed to be constant through out the 
process or made to depend upon gas flow rate itself (as reported in literature [13]). If it is 
assumed to be constant then on substituting Kt in eq. (3.3) we get 

- = -a, • l^-A ■ - + a//J -c (3.4) 

The values of Oi and 7 ] in eq. (3.4) are process dependent and can be obtained by evalua- 
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tion of experimental data. Alternately, if r? is made to depend on gas flow rate itself then 
we get 


dc 

dt 


= -ai ■ 


a dc ^ 


11 a2[-Pi %+QcSf 


(3.5) 


The values of ai and 0-2 in eq. (3.5) are also process dependent and can be obtained 
by evaluation of experimental data. 

In order to determine the best fitting values of ai, o ;2 and rj in eq. (3.4)-(3.5) the 
following approaches have been adopted in the present work for evaluation and/or simpli- 
fication. 


1. Simplification of eq. (3.4) by neglecting the contribution of top jet and bottom gas 
to mass transfer (hereafter referred to as CO-gas evolution model) 

2. Simplification of eq. (3.4) to a non-linear algebraic equation and subsequent evalu- 
ation (hereafter referred to as Non-Linear model) 

3. Simplification of eq. (3.4) to a linear algebraic equation and subsequent evaluation 
(hereafter referred to as Linear model) 

4. Optimization of various models using genetic algorithms (hereafter referred to as 
GA-decarb models) 

5. Purely emperical approach where K' is assumed to be some function of blowing 
time or amount of oxygen blown (hereafter referred to as emperical model) 

6. Statistical models using multiple linear regression (hereafter referred to as Regres- 
sion Models) 

The fundamental basis of each of the above models is described below. 


3.2 CO-Gas evolution model 

If it is assumed that carbon mono oxide gas (CO) evolution rate in the bath is very high 
when compared to the effective top and bottom gas flow rate then the Qe// term in eq. 
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(3.4) can be neglected to give 


dc 

dt 


—a • 



• c 


or 


or 


on integrating, 



= O' • c 



= /3 • ci-i 



dc 

1 



dt 


and assuming that rj is constant, we get 


ct = « + 7i)^ (3.6) 

where (j, = — Y^and optimum value of fj, can be calculated from experimental data by 
using a suitable numerical technique (see chapter 5). 


3.3 Non-Linear model 

Non-Linear model is based on the simplification of decarburization equation (3.4) to a 
non-linear algebraic equation. If it is assumed that stirring energy due to top and bottom 
gas evolution rates is proportional to the stirring energy due to CO evolution* through out 
the process, eq. (3.4) can be simplified to 

'This is also supported by the fact that CO evolution takes place mainly near the impinging spots in the 
hot metal and hence can be assumed to be of similar magnitude to effective stirring energy from top and 
bottom gases . 
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Since ^ < 1 eq. (3.7) can be simplified using binomial expansion to 


dc 

dt 


- - 0:2 • 


« dc 


c 


Upon readjusting the terms 

I [1 - a,o] = -a,c 

Upon integration in the limits Co and Ct we get 




dt 


or 


In - 0:3 • [ct - Co] = -Q !2 • t (3.8) 

Above equation can be further simplified using logarithmic expansion /or— 1 < — l) < 

1 to give 


or 


l-OCz-[Ct- Co] = -012 • t 

Co 


Ct = 


0^2 • Cq 
1 - as 'Co 


• f + Co 


(3.9) 


where 0:2 and 0:3 are constants which can be determined by a suitable numerical procedure 
(see chapter 5). Eq. (3.9) represents a non-linear dependence of q on Co and t and hence 
the name Non-Linear model. 



3.4 Linear model 


Linear model is based on simplification of decarburization equation (3.4) to linear alge- 
braic equation. Equation (3.8) can be rewritten as 


In ct - 0:3 • [ct - Co] = - 0:2 • f + In Co 


Above equation can be further simplified by using logarithmic expansion for — 1 < 
(ct - 1 ) < 1 and -1 < (co - 1 ) < 1 to 


Ct — Otz- [Ct — Co] = — Q ;2 • f -I- Co 


or 


ct = -p-^-f + Co (3.10) 

[1 - 0:2] 

where 0:2 and az are constants which can be determined by a suitable numerical procedure 
( see chapter 5). Eq. (3.10) represents a linear dependence of ct on Cq and t and hence the 
name Linear model. 


3.5 GA-decarb models 

These models use genetic algorithms [18] to dynamically evaluate the optimum mass 
transfer coefficient at each time step. Two methods can be followed to dynamically opti- 
mize eq. (3.4). 

Method 1 

In this method, rj is assumed to be constant through out the process. Eq. (3.4) can be 
written as 


Ac = At • O'! • 


^ dc ^ 


c 


(3.11) 
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For small time steps of At, the parameters to be optimized are ai and ri and are determined 
by using GA (as explained in chapter 5). Carbon at the present time step, Ct, is calculated 
as 


ct = Co - Ac 

where Co is the carbon content at previous time step. 

Optimization is done successively for time steps ( say At = O.Olsec) until t < tmax 
where tmax is the total time of second blow period. 

Method 2 

Since the gas flow rate during the process is changing substantially, rj can be made to 
depend on the gas flow rate itself. The governing equation for optimization becomes 

/ dc 

Ac = Af • Q!X • f --/5i + Qeffj ■ c (3.12) 

For small time steps of At, the parameters to be optimized are ai and a 2 . These are 
determined by using GA (as explained in chapter 5). 


3.6 Emperical model 

Since gas evolution rate slows down with time and time of second blow is short, capacity 
mass transfer coefficient K' in eq. (3.3) can be made to depend inversely on the time of 
second blow period t as 



or 
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where, C is a proportionality constant. Upon substitution of k in (3.3) , we get 

dc _ c 
dt t 

On integration, 


dc 
J Co C 



or 


or 






(3.13) 


Optimum value of ^ can be calculated from experimental data (see chapter 5). 


3.7 Regression models 

The models discussed above can be further subjected to multiple linear regression (MLR) 
to incorporate the effect of other operational variables on the end point carbon prediction. 
For example, for the case of Linear model, the equation for end point carbon prediction 
can be written as 


n—l 

4 = Mo + Ml Co + X) (3-14) 

i=2 

where, is the predictedend point carbon and Co is the carbon value at sublance from 
actual model and Mi ’s are coefficients (to be determined by MLR) of operational variables 
Xi ’s such as basicity, orel^, dolol, hlans2, ore2, dolo2, slag2 etc. It may be noted that 
^These variables are explained earlier. 
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regression model is nothing but a special case of Linear model discussed earlier (Section 
3.4). 

Regression models can also be developed for prediction of temperature and dissolved 
oxygen by incorporating all the variables which affect the process 

n— 1 

r - fio + ‘ (3.15) 

1=1 

where r = T for Temperature prediction models 
r = [O] for oxygen prediction models. 

fii ’s are coefficients (to be determined by MLR) of operational variables ’s such as 
basicity, Orel, dolol, Ti , hlans2, ore2, dolo2, slag2 etc. For details of variables see List 
of Symbols. 


3.8 Models for phosphorus prediction 

At the end of blow in EOF steel making only a quasi static (or pseudo thermodynamic) 
equilibrium is approached between slag and metal with a separating interface. Since the 
slag only approximates to equilibrium with metal and gas[CO, CO 2 ] at the end of blow, 
the application of thermodynamics for direct estimation of phosphorus distribution at end 
point has not been straight forward. The models developed to predict endpoint phospho- 
rus distribution have been generally based on regression analysis of actual plant data by 
assuming a pseudo thermodynamic equilibrium at the end of process. In the present work, 
however, a kinetic model has also been tried. 

3.8.1 Thermodynamics of phosphorus distribution 

The molecular theory of slag stipulates the formation of tricalcium phosphate in lime rich 
slags 

[P]+5/2[FeO]+3/2CaO = l/2{Ca3 (P04)2}+5/2(Fe) 


38 



(3.16) 


K - (PO4) 

~ „ 1.5 . „2 5 . 

°CaO ®FeO “f" 

This process can be vied as a multi step process involving several stages 

i) [Pm] Pg (gas), at slag/metal interface 

ii) P(p) + 0(p) -4- P205(g) at slag/metal interface 
ii) P205(5) + {CaO) —¥ {CaO .P^O^') 

Step (ii) in the gas-slag (transitory intermediate) reaction introduced by Wagner[17 ] 

as 

(P2), + 5/2 (O2), + 3 (O) -» 2 (PO4)-’ 

and is the basis of the term phosphate capacity, a measure of dephosphorising capability 
of slag under a given oxygen potential and is related to oxygen activity in slag by 

g} = C {PO<)-= • 

where, C(P04)~® is the phosphate activity which is strongly influenced by temperature 
and slag composition. The conversion of P to P2O5 takes place at slag metal interface and 
is related to activity of P in metal and other active parameters of slag as implied above. 
The rate determining stage of the reaction is activity of [P] at the slag metal interface and 
is decided by the rate of transport of phosphoras to the slag metal boundary for oxidation. 
If manganese is present in appreciable concentrations then dephosphorization process is 
hindered. 

Balajiva et al [14] proposed dephosphorisation reaction as 

2P[metal]+5Feo[slag] = P205[slag]+Fe[metal] 

for which the apparent equilibrium constant can be written as 

_ Wt%P 205 

~ [wt%P]'^[wt%FeO]^ 

Regression analysis of the above equation on data obtained during laboratory experiments 
yielded 

log K = 10.78%CaO-0.089T-6.425 where T is Deg C. 

In the above study the three important variables identified for dephosphorisation were 
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(i) FeO content of slag which represents oxidation potential of slag, 

(ii) CaO content of slag which is directly related to basicity and 

(iii) Turn-down temperature [T], deg C and a negative coefficient for temperature im- 
plies poorer dephosphorization as temperature rises. 

The simultaneous effect of temperature, lime and FeO content of slag on phosphoms 
distribution has been described by Healy [14] by using ionic theory of slag 


log wCf = + 0.08(%CaO) + 2.5 log(%Fe) - 16.0 

[/of\ 1 


(3.17) 


where [%P]=mass % of phosphoms in metal 
(%P)=mass % of phosphorus in slag 
(%Fe)=mass %Fe in slag 
(%C7oO)=mass %CaO of slag 
r=temperature in 

Healy’s equation is widely used in steel industry by adapting it to the plant data. 

Since temperature at tap varies in a small range, Deo et al [15] proposed a dephospho- 
rization model in terms of basicity only (on the basis of plant data), 
log( %P)/[P] = 0.43-Basicity 0.55 
where, basicity = %Ca0/%Si02 


3.8.2 Phosphorus prediction models developed in the present work 

Phosphorus prediction is done using 

1 . Kinetic model of mass transfer 

2. Regression model 

3.8.2.1 Kinetic model of mass transfer 

Similar to the decarburization equation (3.1), flux of phosphoms in metal is given by 

f = (3.18) 
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or 

| = -A'V-P-) 

where kp is mass transfer coefficient for phosphorus, A is area of slag-metal interface, V 
is volume of bulk metal (including droplets), p* is concentration of phosphorus in the bulk 
and p* is the concentration of phosphoms at the interface and K' is capacity mass transfer 
coefficient. If it is assumed, as also done in the previous works[ ], that mass transfer 
coefficient is nearly constant through out the process then eq. (3.18) can be integrated to 
give 




Po-Pe 


Vpo-Pe 


(3.19) 


where t is time in seconds, pt is end point phosphorus, Po is sublance phosphoms and Pe 
is the equilibrium concentration of phosphoms. Eq. (3.19) can be simplified to 


In ^* (3.20) 

Po-Pe 

However, K' is dependent upon operational parameters such as basicity, hot ratio, 
Orel, dolol, slagl, hlans2, ore2, dolo2, slag2 etc. The effect of other variables affecting 
phosphoms distribution can be incorporated by MLR and hence above equation can be 
written as 


- In ^ = fio+ E l^^ • (3-21) 

Po - Pe ^ 

where pt ’s are coefficients of operational variables Xi such as basicity, hot ratio, orel, 
dolol, slagl, hlans2, ore2, dolo2, slag2 etc. 

3.8.2.2 Regression model 

This model is similar to the one developed for carbon (eq. (3.14)). Here end point phos- 
phoms is predicted as a function of operational variables. 

n— 1 

Pt — Po A piPo + E (3.22) 

1=2 
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where, p[ is the end point phosphorus value predicted from MLR, Po is sublance phospho- 
rus and Pi ’s are coefficients (to be determined by MLR) of operational variables Xj such 
as basicity, hot ratio, orel, dolol, slagl, hlans2, ore2, dolo2, slag2 etc. It may be noted 
that regression model is only a special case of eq. (3.21). 


3.9 Models for manganese prediction 

The FeO and MnO exchange reaction involving the oxidation of manganese in steel can 
be formulated as 


{FeO) + [Mn] = {MnO) -f [Fe] 

The equilibrium relation can therefore described in terms of the mass concentration of 
oxides 

{%MnO) 

^FeMn (%FeO)[%Mn] 

where the equilibrium relation K" depends on temperature and slag composition. 

3.9.1 Thermodynamics of manganese distribution 

TurkdogaCn[16] proposed the relationship for manganese partition as 
(MnO%)/[%Mn] = 6(%Feo)*Si02/CaO, 

where combined effect of slag FeO and basicity of Mn partition have been considered. 
Two new relations for estimating Mn partition have also been proposed in which 

(i) hot metal Si is considered, since it remarkably affects basicity (slag route) and FeO 
in slag. 

(ii) residual oxygen activity of the steel, which is related to the end point carbon, 
has been introduced. This residual oxygen does not participate in any of the oxidation 
reactions in the melt, but it may activate the manganese atoms effectively such that the 
kinetics of Mn + FeO = MnO + Fe and subsequent removal of MnO to slag phase is 
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enhanced at the slag/metal interface. In other words, higher the oxygen activity of the 
bath, larger will be Mn partition to slag i.e., Mn recovery will be less. 

Deo et al[15] proposed a model for manganese partition as 
(Mn)/[Mn] = 2.1864(%Fe203+%FeO)-23.728 
and since at the quasi-static stage (FeaOs+FeOlt can be expressed in terms of FeOt, the 
manganese partition is expressed as 
(Mn)/[Mn] = 2. 1 864%FeO-23.728, 

where the only variable is slag FeO. This equation is, however, applicable only to a 
narrow and specific range of operation. 

3.9.2 Models developed for manganese prediction in the present work 

Manganese prediction was done using 

1. Kinetic model of mass transfer 

2. Regression model 

3.9.2.1 Kinetic model of mass transfer 

Similar to decarburization equation (3.1), flux of manganese in metal is given by 

^ = -kMu^{Mn^ - Mn^) (3.23) 

or 


^ - Mn‘) 

at 

where kMn is mass transfer coefficient for manganese, A is area of slag-metal interface, V 
is volume of bulk metal (including droplets), Mn!’ is concentration of manganese in the 
bulk and Mn* is the concentration of manganese at the interface and K' is capacity mass 
transfer coefficient. If it is assumed , as done in the previous works[17], that K' is constant 
through out the process then eq. (3.23) can be integrated to give 
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(3.24) 


In 


Mrit — Mue 




Mn„ 


Mrio - Mrie V Mrio - Mn^ 

where Mrit is end point manganese, Muo is sublance manganese and Mn^ is the equilib- 
rium concentration of manganese. Eq. (3.24) can be simplified to 


, Mrit — Mrie 
^ Muo - Mne 


-K't 


(3.25) 


However, K' is dependent upon operational parameters such as basicity, hot ratio, 
Orel, dolol, slagl, hlans2, ore2, dolo2, slag2 etc. Other variables affecting manganese 
distribution can incorporated by MLR and hence above equation can be written as 


-In 


Mrit — Mrie 

Mrio — Mrif 


n-l 

— "b ^ ^ l^i ■ 
1=1 


(3.26) 


where /ij ’s are coefficients of operational variables such as basicity, hot ratio, orel, 
dolol, slagl, hlans2, ore2, dolo2, slag2 etc. 


3.9.2.2 Regression model 

This model is similar to the one developed for carbon (eq. (3.14)). Here end point man- 
ganese is predicted as a function of operational variables. 

n-l 

Mti[ = /Xo + fiiMuo + (1% • Xi (3.27) 

2=2 

where, Mn't is the end point manganese value predicted from MLR, Mrio is the sublance 
manganese and jXi ’s are coefficients of operational variables such as basicity, hot ratio, 
Orel, dolol, slagl, hlans2, ore2, dolo2, slag2 etc. It may be noted that regression model is 
only a special case of eq. (3.26). 
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Chapter 4 

Experimental Data 


Experimental data was collected for three converters of 300, 200 and lOOton capacity 
from literature and is classified as plantl [1 8], plant2[19] and plant3[20], respectively. The 
following section explains the data and classification. 


4.1 Plantl 

Table 3.1 shows a sample set of heats for a 300 ton converter consisting of parameters 
basicity, hotratio, hot ratio*initial phosphoras (p), slag volume, ore added during first 
blow, slag added during first blow, dolomite added during first blow, sublance temperature 
(To), sublance carbon (Cq), oxygen blown during second blow (022), slag added during 
second blow (slg2), dolomite added during secon blow (dolo2), lance height during second 
blow (hlans2), end point carbon (ct), sublance phosphoms (po), end point temperature (T), 
end point manganese (Mn*), end point phosphorus (po), sublance manganese (Muq) and 
end point dissolved oxygen ([O]) used in tuning the models. Oxygen gas flow rate is 750 

3 

The parameters calculated from actual heat are defined below, for a typical heat with 
hot metal weight (HMwt) = 295400 kg, scrap weight = 56700 kg, Silicon content of hot 
metal (%Si) = 0.359%, lime added during first blow (lime 1) = 1 1956 kg, dolomite added 
in the first blow (dolol)=413 kg and phosphoms content of hot metal (%p)=0.056 %, 
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= 2272.47 Kg. 


c /n xjT^f 4. 

Si02 = HMwt * — — ★ — 
100 28 

0-359 60 

= 295400 ■*■ ~ -k — 
100 28 


Slag volume = Si02 +lime l|slgl = 2272.47+1 1956+2980 = 17208.47 Kg 


hot ratio = 


HMwt 295400 


scrapwt 56700 


= 5.2098 


basicity = 


limel + dolol 


11956 + 413 
2272.47 


5.443 


The actual data is divided into 3 case studies. 


4.1.1 Case study 1 

Number of heats in this case study are 273. For all the modes developed on these 273 
heats, it is referred to as case study 1(273 heats ). Two subsets of this data are created 
by subjecting it to the overall restriction that end point carbon should be in the range of 
0.035 %-0.075% and whether ore is added or dolo2 is added during second blow. 

4.1.2 Case study 2 

End point carbon is in the range of 0.035-0.075% and ore 2 is not added in the heats of 
this case study. Number of heats in this case study are 156. For the models tuned on these 
156 heats, it is referred to as case study 2 ( 156 heats ). 

In 33 heats Ore2 is added in place of dolo 2 in the heats. Since number of heats in this 
group are only 33 no models are tested on this set. 
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Heat No Slag volume hotratio *p* Po T Pt Mug Mut [O] hlans2 
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M9250 25058.34 0.528 0.023 1702 0.008 0.290 0.140 595 180 0.046 

M9319 1 7208.47 0.292 0.027 1662.000 0.008 0.326 0.174 413 187 0.050 

Mean 17886.256 0.309 0.022 1 665 0.009 0.316 0.176 426 194 0.059 

a 1857.449 0.041 0.005 14 0.002 0.029 0.026 89 22 0.011 



4.2 Plant 2 


Table 3.2 shows a sample set of heats for a 200 ton converter consisting of Si02, 022, 
scrap added during first blow, sublance temperature (To), end point temperature (T), sub- 
lance carbon (Co), flux added, waste oxide briquettes (WOBS), chill scrap, end point car- 
bon {ct) and end point dissolved osygen ([O]). St02 is calculated as explained earlier for 
plant 1. Oxygen gas flow rate is 600^. Heats are selected with end point carbon lying 
in the range of 0.02 to 0.075% only. There is a total of 321 heats. Parameter SiO^ is 
calculated as already shown for plant 1. 


4.3 Plant 3 

Table 3.3 shows a sample set of heats for a 100 ton converter consisting of parameters 
basicity, slag volume, sublance temperature (To), sublance carbon (Co), oxygen blown 
during second blow (022), raw dolomite added during second blow (raw dolo2), burnt 
dolomite added during second blow (burnt dolo2), lance height during second blow (lance 
level), end point carbon (ct), sublance phosphorus (po), end point temperature (T), end 
point manganese (Mn*), end point phosphorus (po). sublance manganese (Mn^) and end 
point dissolved oxygen ([O]) used in tuning the models. End point carbon range for these 
heats are 0.03 to 0.09% only. Oxygen gas flow rate is 340^. Heats are selected with 
end point carbon lying in the range of 0.03 to 0.07% only. There is a total of 297 heats. 
Method of calculation of parameters such as basicity, slag volume is same as for plant 1 . 

Comparison of mean and cr values for the variables of Plant 1, Plant 2 and Plant 3 is 
shown in Table 4.4. 
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Table 4.2 Sample data set for Plant 2 
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Table 4.4; Comparison of Mean and <r for operational variables of Plant 1, Plant 2 and 
Plant 3. 


Operational 

Variavles 

Plant 1 

Mean a 

Plant 2 

Mean a 

Plant 3 

Mean a 

Basicity 

5.054 

0.831 

- 

- 

7193.807 

1461.411 

hot ratio 

5.305 

0.979 

- 

- 

- 

- 

Orel (kg) 

4782 

1627 

- 

- 

- 

- 

slgl (kg) 

3021 

223 

- 


- 

- 

dolol (kg) 

819 

626 

- 

- 

- 

- 

To {°C) 

1603 

17 

1584 

25 

1629 

21 

Co (mass%) 

0.412 

0.082 

0.362 

0.182 

0.385 

0.134 

022 

2149 

304 

1402 

361 

557 

156 

Ore 2 (Kg) 

116 

283 

- 

- 

- 

- 

slg2 (kg) 

91 

372 

- 

■■ 

- 

- 

dolo2 (kg) 

540 

553 


- 

- 

- 

Slag volume (kg) 

17886.256 

1857.449 

- 

- 

3137.791 

460.325 

hotratio * p (kg) 

0.309 

0.041 


- 

- 

- 

Pomass% 

0.022 

0.005 

- 

- 

0.015 

0.004 

T("C) 

1665 

14 

1661 

16 

1689 

19 

Pt mass% 

0.009 

0.002 


- 

0.008 

0.002 

Muo mass% 

0.316 

0.029 

- 

- 

0.202 

0.0041 

Mrit mass% 

0.176 

0.026 

- 

- 

0.144 

0.026 

[0] ppm 

426 

89 ^ 

893 

296 

602 

161 

hlans2 cm 

194 

22 

- 

- 

159 

3 

ct mass% 

0.059 

0.011 

0.044 

0.019 

0.051 

0.010 

Si02 

- 

- 

2083.849 

575.331 

- 

- 

scrap 

- 

- 

35407 

9492 

- 

- 

chill scrap (Kg) 

- 

- 

335 

799 

- 

- 

wobs (Kg) 

- 

- 

4774 

2001 

- 

- 

flux (Kg) 

- 

- 

9878 

1954 

- 

- 

raw dolo2 (Kg) 

- 

- 

- 

- 

247 

310 

Burnt dolo2 (Kg) 

- 

- 

- 

- 

12 

106 

Lance life 

- 

- 

- 

- 

124 

70 
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Chapter 5 


Results and Discussion for End point 
Carbon, Temperature, Dissolved 
oxygen. Phosphorus and Manganese 
Models 

5.1 Introduction 

This chapter discusses the results obtained for various end point carbon, temperature, 
phosphorus, manganese and dissolved oxygen models presented in chapter 3. The ex- 
perimental data for Plant 1, Plant2 and Plant3, on which models heve been tested and 
tuned, are already presented in chapter 4, 

5.2 Results of Plant 1 

Two separate data sets have been considered for Plant 1 . As described earlier, for Plant 
1 the set of 273 heats is designated as Case Study 1. A sub set consisting of 156 heats, 
where only those heats in which raw dolomite is added as a coolant during second blow, 
is called Case Study 2. The results of all models for Plant 1 are summarized in Table 5.1 
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and results obtained with each model is are presented below, individually. 


5.2.1 End point carbon prediction models for Plant 1 

Fundamental aspects of development of end point carbon prediction models are discussed 
in chapter 3. 

5.2.1.1 CO-Gas Evolution model (model 5.1) 

CO-Gas evolution model is given by eq. (3.6) 

ct = {c^ + rt)^ (5.1) 

where /i = model is developed for Case Study 2 (156 heats) of Plant 1. 

By using Gauss-Newton method the sum of square of errors (SSE) (]2(actual— predicted)^) 
was obtained for different values of ju in the range -1 .5 to -0.33 (i.e. 77 in the range of 0.25- 
0.6). The optimum value of /x calculated from the graph of sum of square of errors (SSE) 
vs n (Fig. 5.1, corresponding to minimum in SSE is) -0.41. The corresponding r) value 
can be back calculated as 0.291 . The 7 value, corrsponding to minimum SSE, is 0.010292. 

Fig 5.2 shows the graph of actual vs predicted carbon. From the statistics of regression 
performed on actual vs predicted carbon ( see Table 5.1, model 5 . 1 ), it can be observed that 
a value is ±0.007 (0.0066)%C and correlation coefficient R is 0.59 (0.5852). Therefore 
model 5.1 is able to explain only 34% of variation. For a typical heat with 
Co = 0.517%, 022 = 2462, t = 197s, we obtain by substituting in eq. (5.1) 
ct = (0.517-o-^^± 0.01029 * 197) ^ 

= 0.052. 

Actual end point carbon value for this particular heat is 0.050. 


5.2.1.2 Non-Linear model (model 5.2) 

Nonlinear model is given by eq. (3.9) 


Ct = 


0:2 • Cq 
1 — 0!3 • Co 


• t + Cq 


(5.2) 
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Table 5.1 Summary of results for carbon prediction obtained for Plant 1 (contd... 
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Table 5.1 Summary of results for carbon prediction obtained for Plant 1 (contd... 
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(reg) Z?Z 2 l^i-Xi dolol.To (0.060) (0.050) c=0.0261 study 2 

hlans2 R=0.74 

Co, t (7=0.004 

F=1 82.483 



Sum of square of errors 



- 1.6 - 1.4 - 1.2 -1 - 0.8 - 0.6 - 0.4 - 0.2 

mu 

Figure 5.1; mu (n) vs sum of square of errors for CO-Gas evolution model (model 5.1); 
minimum SSE occurs at /x = -0.41. 
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Predicted carbon (mass%) 



Actual carbon (mass%) 


Figure 5.2: Data from Plant 1, Case Study 2 (156 heats): Actual vs predicted carbon for 
CO-Gas evolution model (model 5.1), best fit line is y=0.5694x+0.0251. Statistics are 
n=156, R=0.59, std. error=0.007. 
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where q ;2 and 0-3 are constants to be tuned. This model is tuned for Case Study 2 (156 
heats) of Plant 1 . Constants 0:2 and 03 , again determined by applying Gauss-Newton 
method on 156 heats, are found to be 0.0075 and -1 .2654 , respectively. From the statistics 
of regression performed on actual versus predicted carbon it can be observed that <r=0.024 
(0.02449 )%C and R=0.66 (0.6557) (see Table 5.1, model 5 . 2 ). Thus, this model is able to 
explain 43% of variation. For a typical heat with cq = 0.517%, t = 197 s we get 

_ -0.0075*0.517*197 | ni=;i7 
^ ~ 1+1.2654*0 517 ' 

= 0.055. 

Actual end point carbon value for this particular heat is 0.050. 

5.2.1.3 Linear model (model 5.3) 

Linear model is given by eq. (3.10) 

Q!2 

Ct — “r; 7 ’t + Co 

[1 - as] 

where 0:2 and 0:3 are constants. The above equation can be simplified into the form 

Cf = —0 • f ■+• Co (5.3) 

where the constant 0 is determined by performing linear regression on eq. (5 .3) as 0.002062, 
When this model is tuned for 156 datasets of Plant l,it can be observed from the statis- 
tics of regression performed on actual versus predicted carbon that cr=0.04 (0.0397)%C 
and R=0.4935 (0.49) (see Table 5.1, model 5.3). Thus, this model is able to explain only 
24.35% of variation. For a typical heat with t=197 s, Cq = 0.517%, we get from eq. (5.3) 
ct =- 0.002062 * 197 + 0.517 
= 0 . 11 . 

Actual end point carbon value for this particular heat is 0.050. 

5.2.1.4 GA-decarb models (model 5.4) 

GA-decarb models have been discussed earlier in section 3.5. In these models, genetic 
algorithm (GA) technique is used to optimize the parameters. After trail and error runs, 


62 



the GA parameters used for optimization of these models are: micro GA(Yes) [?], Cross 
over probability (0.9), Creep probability (0.004). Mutation probability (0.00001), Number 
of generations (300), Population size (10), Chromosome size (10). GA code developed by 
David L. Carroll of University of niinois[?] was used for optimizing the models. 

Method 1 

The governing equation for optimization is given by eq. (3.1 1) 

= + -c (5.4) 

The parameters to be optimized are oii and 77 . Optimization is done successively for small 
time steps ( say At = O.Olsec) until t < tmax where tmax is the total time of second blow 
period. Range of ai is 0.005-0.007 and range of 77 is 0.20-0.65. Optimum values of ai and 
77 are selected coixesponding to minimum sum of square of eirors (SSE) between actual 
and predicted carbon. 

Case Study 1 (273 heats) 

Eq. (5.4) was optimized for 273 data sets. Optimized values of ai and 77 corresponding to 
minimum SSE are 0.006771 and 0.2247, respectively. Table 5.1, model 5.4 (Case Study 
1) shows the statistics of regression performed on actual vs predicted carbon. It can be 
observed from regression statistics (Table 5.1, model 5.4) that a=0.008 (0.00806)% and 
R=0.55 (0.5552), thus explaining 30.8% of variation. Fig. 5.3 shows graph of actual vs 
predicted carbon. The significance of variation of mass transfer coefficient {K'), where 
K' = ai ■ + Qeff) . is discussed later. 

Case Study 2 (156 heats) 

Eq. (5.4) was optimized for 156 data sets. Optimized values of ai and 77 corresponding to 
minimum SSE are 0.006501 and 0.2477. The statistics of regression performed on actual 
versus predicted carbon Table 5.1, model 5.4 (Case Study 2) shows that <7=0.007 (0.0071) 
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Predicted carbon (mass %) 


0.07 


Actual carbon (mass%) 

Figure 5.3: Data from Plant 1, Case Study 1(273 heats): Actual vs predicted for model 
4 (Case Study 1), best fit line is y=0.4987x+0.0302. Statistics are n=273, R=0.55, std. 
error=0.008. 



%C and R— 0.60 (0.6015), thus explaining 36.2% of variation. Fig. 5.4 shows graph of 
actual vs predicted carbon. 

Method 2 

In this method, rj is made proportional to gas flow rate, because gas flow rate changes 
substantially with time. The governing equation for optimization is given by eq. (3.12) 
i.e.. 


Ac = Af • O'! 


dc 

dt 


Pi + Qeff 


“2 


(5.5) 


The parameters to be optimized are ai and 02 . Optimization is done successively for timp. 
steps ( say At = O.Olsec) until t < t^ax where tmax is the total time of second blow 
period. Range of cti is 0.01-0.016 and range of 0:2 is 0.0005-0.0013. Optimum values of 
ai and 0:2 are selected corresponding to minimum sum of square of errors (SSE) between 
actual and predicted carbon for all heats. 


Case Study 1 (273 heats) 

Eq. (5.5) was optimized for 273 data sets. Optimized values of ai and 0:2 correspond- 
ing to minimum SSE are 0.01 109 and 0.0006283, respectively.The statistics of regression 
performed on actual vs predicted carbon, Table 5.1, model 5.4 (method 2, Case Study 1) 
shows that cr=0.010 (0.00947)% and R=0.62 (0.6180), thus explaining 38.2% of variation. 
Fig. 5.5 shows graph of actual vs predicted carbon. 


Case Study 2 (156 heats) 

Eq. (5.5) was optimized for 156 data sets. Optimized values of ai and a 2 corresponding 
to minimum SSE are 0.01 105 and 0.0007495. From the statistics of regression it can 
be observed that c7-=0.008 (0.00826)%C and R=0.67 (0.6693) (see Table 5.1 model 5.4 
(method 2, Case Study 2)), thus explaining 44.8% of variation. Fig. 5.6 shows the graph 
of actual vs predicted carbon. 
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Predicted carbon (wt%) 



Actual carbon (wt%) 

Figure 5.4: Data from Plant 1, Case Study 2 (156 heats): Actual vs predicted graph for 
model 5.4 (method 1, Case Study 2), best fit line is 0.6292x+0.0216. Statistics are n=156, 
R=0.60, std. error=0.007. 
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Predicted carbon (mass %) 



Actual carbon (mass%) 

Figure 5.5: Data from Plant 1, Case Study 1(273 heats): Actual vs predicted carbon for 
model 5.4 (method 2, Case Study 1 ). Best fit line is y=0.6898x+0.0182. Statistics are 
n=273, R=0.62, std. error=0.010. 
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Figure 5.6: Data from Plant 1, Case Study 2(156 heats): Actual vs predicted graph for 
model 5.4 (method 2, Case Study 2), best fit line is y=0.8781x+0.0076. Statistics are 
n=156, R=0.67, std. error=0.008. 




5.2.1.5 Empirical Model (model 5.5) 

Empirical model is given by eq. (3.13) i.e., 



(5.6) 


Optimum value of C is calculated from the graph of sum of square of errors (SSE) vs 
C- ^ typical graph of SSE versus ^ is shown in Fig. 5.7. From the graph it can be observed 
that optimum value of C corresponding to minimum SSE is 0.382. From the statistics of 
regression it can be observed that cr=0.009 (0.0089)% and R=0.37 (0.3663) (see Table 5.1 
model 5.5), thus explaining only 13.4% of variation. Fig 5.8 shows the actual vs predicted 
carbon. For a typical heat with Cq = 0.517% and t=197 s, we get 


=0.068. 

Actual end point carbon for this particular heat is 0.050. 


5.2.1.6 Regression models 

As an attempt to improve upon the prediction of models developed so far, the effect of 
other variables is incorporated in regression models. Regression models as discussed in 
section 3.7 are given by eq. (3.14) i.e., 

n-l 

Cj = fJ-o "b fJ'lCo + (5.7) 

1=2 

where, cj is the end point carbon and Co is the spblance carbon value, /ij ’s are coefficients 
of operational variables X^. These models are developed using multiple linear regression 
(MLR) technique. The units of the operational variables are as shown in Table 4.1 


Model 5.1 (reg) 

Model 5.1 (reg) is regression model for CO-Gas evolution model (model 5.1). Variables 
accepted and corresponding coefficients in this model are Orel (-0.0000009), dolol(- 
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Predicted carbon (mass %) 


0.075 

0.07 

0.065 

0.06 

0.055 

0.05 

0.045 

0.04 

0.04 0.045 0.05 0.055 0.06 0.065 0.07 0.075 

Actual carbon (mass %) 

Figure 5.8: Data from Plant 1, Case Study 2(156 heats): Actual vs predicted carbon for 
model 5.5. Best fit line is y=0.412x+0.034. Statistics are n=156, R=0.37, std.erroi^.008. 
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0.0000037), Cf (0.5438315) and constant(0.0317). From the statistis of regression it can 
be observed that cr=0.004 (0.00417)%C and R=0.66 (0.6571) (Table 5.1, model 5.1 (reg)), 
thus explaining 43.2% of variation. Fig. 5.9 shows the graph of actual versus predicted 
carbon. It can be observed that R increased from 0.59 in model 5.1 to 0.66 in model 5.1 
(reg) i.e., an improvement of 12.3%. The value of o also came down to 0.004% in model 
5.1 (reg) from 0.007% in model 5.1, a decrease of 36.8% in standard error. For a typical 
heat with ore 1 = 3994, dolol = 413 and q = 0.052 (see model 5.1), 
c[ = 0.0317 + 0.5438315(0.052) - 0.0000009(3934) - 0.0000037(413) 

= 0.055. 

The actual carbon for this heat is 0.050. 

Model 5.2 (reg) 

Model 5.2 (reg) is regression performed on Non-Linear model (model 5.2). The respec- 
tive coefficeints of operational variables by MLR acceped are Orel (-0.0000007), dolol(- 
0.0000026), Ct (0.1644613) and iXo (0.0521). From the statistis of regression (Table 5.1 
(model 5.2 (reg))) it can be observed that (j=0.004 (0.00396)%C and R=0.69 (0.6847), 
thus explaining 46.9% of variation. Fig. 5.10 shows the graph of actual versus predicted 
carbon. It can be observed that R increased from 0.66 in model 5.2 to 0.69 in model 5.2 
(reg) i.e., an improvement of 4.5%. a also came down to 0.004% in model 5.2 (reg) from 
0.024 in model 5.2 i.e., a reduction of 83.8% in standard error. For a typical heat withorel 
= 3994, dolol =413 and Ct = 0.055 (see model 5.2), 

<4,^0.0521 + 0.1644613(0.055) - 0.0000007(3994) - 0.0000026(413) 

= 0.057. 

The actual carbon for this heat is 0.050. 

Model 5.3 (reg) 

Model 5.3 (reg) is regression model for Linear model (model 5.3). Operational parameters 
accepted in this model are Orel (-0.0000007), dolol (-0.0000022), 022 (-0.0000227), 
hlans2 (0.0000728), Ct (0.0843022) and jU^ (0.0623). From the statistis of regression (Table 
5.1, model 5.3 (reg)) it can be observed that cr=0.004 (0.0041 8)%C and R=0.76 (0.7617), 


72 



Predicted carbon (wt%) 


0.075 


0.07 


0.065 


0.06 


0.055 


0.05 


0.045 


0.04 


0.035 

0.035 0.04 0.045 0.05 0.055 0.06 0.065 0.07 0.075 

Actual carbon (wt%) 

Figure 5.9: Data from Plant 1, Case Study 2(156 heats): Actual vs predicted carbon for 
model 5.1 (reg), best fit line is y=0.4294x+0.033. Statistics are n=156, R=0.66, std. er- 
ror=0.004. 
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Figure 5.10: Data from Plant 1, Case Study 2(156 heats): Actual vs predicted carbon 
for model 5.2 (reg), best fit line is y=0.4980x+0.0315. Statistics are n=156, R=0.69, std. 
error=0.004. 
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thus explaining 58% of variation. Fig. 5.11 shows the graph of actual vs predicted carbon. 
It can be observed that R increased from 0.5 in model 5.3 to 0.76 in model 5.3 (reg) i.e., 
an improvement of 54.3%. The value of a also came down to 0.004% in model 5.3 (reg) 
from 0.04% in model 5.3 i.e., a reduction of 100% in standard error. For a typical heat 
with Orel = 3394, dolol = 413, 022 =2464, hlans2 = 187, we get 

c'j = 0.0623 + 0.0843022 (0.517) - 0.0000007(3994) - 0.0000022(413) - 
0.0000227(2464) + 0.0000728(187) 

= 0.059. 

The actual carbon for this heat is 0.050. 

Model 5.3 (a) (reg) 

In this model sublance phosphoms and manganese are also considered as parameters in 
MLR. From statistics of regression (Table 5.1, model 5.3 (a) (reg)), it can be observed 
that R=0.82 and (j=0.004. Though there is no improvement in a value from model 5.3 
(reg) to model 5.3 (a) (reg), there is 7.7% improvement in R value. Thus if sublance 
phosphorus and manganese are available, this model gives better results. Fig. 5.12 shows 
graph of actual vs predicted carbon. Variables used in this model are orel(-0.0{X)0008219), 
ro(0.000 1869 125), Co(0. 1146421219), Muo (0.0843178237), po(-l. 0100703893), 022(- 
0.0000249837) and constant (-0.2346). For a typical heat with orel=3994. To =1607, Co - 
0.326, Mno=0.517, Po=0.027, 022=2464 

c't=-0.2346-0.0000008219(3994)+0.0001869125(1607)+0.1146421219(0.326)+ 

0.0843178237(0.517)-1.0100703893(0.027)-0.0000249837(2464) 

=0.055. 

Actual end point carbon for this heat is 0.050. 

Model 5.4 (reg, method 1, Case Study 1) 

Model 5 .4(reg, method 1 , Case Study 1 ) is regression model for model 5 .4 (method 1 , Case 
Study 1). The operational variables effecting the end point carbon content and statistics of 
regression are presented in Table 5.1, Model 5.4(reg, method 1, Case Study 1). Variables 
used in this model are hot ratio (0.0027146), Orel (-0.0000014), ct(0.6304038), Dolol(- 
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Figure 5.11: Data from Plant 1, Case Study 2(156 heats): Actual vs predicted carbon 
for model 5.3 (reg), best fit line is y=0.5796x+0.0245. Statistics are n=156, R=0.76, std. 
error=0.004. 
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Figure 5.12: Data from Plant 1, Case Study 2(156 heats): Actual vs predicted carbon for 
model 5.3 (a) (reg), best fit line is y=0.6708x+0.0192. Statistics are n=156, R=0.82, std. 
error=0.004. 
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0.0000029), 0re2(-0.0000063), To(-0.0000870) and Constant(0.1542). Fig. 5.13 shows 
the graph of actual vs predicted carbon. It can be observed that R increased from 0.55 in 
model 5.4 (reg) to 0.67 in model 5.4 i.e., an improvement of 20%. The value of a also 
came down to 0.005 in model 5.4 (reg, method 1, Case Study 1) from 0.008 in model 5.4 
(method 1, Case Study 1) i.e., a reduction of 33.6%. 

Model 5.4 (reg, method 1, Case Study 2) 

Model 5.4 (reg, method 1, Case Study 2) is regression model for model 5.4 (method 1, 
Case Study 2). Variables accepted in this model are Orel (-0.0000009), Dolol (-0.0000036), 
Cf(0.5167197) and Constant(0.0334) . Statistics of regression are presented in Table 5.1, 
model 5.4 (reg, method 1, Case Study 2). Fig. 5.14 shows the graph of actual vs predicted 
carbon. It can be observed that R increased from 0.60 in model 5.4 (method 1, Case Study 
2) to 0.66 in model 5.4 (reg, method 1, Case Study 2) i.e., an improvement of 9.7%. The 
value of a also came down to 0.004 in model 5.4 ( reg, metiiod l,Case Study 2) from 0.007 
in model 5.4 (method l,Case Study 2) i.e., a reduction of 41% . 

Model 5.4 (reg, method 2, Case Study 1) 

Model 5.4 (reg, method 2, Case Study 1) is regression model for model 5.4 (method 
2, Case Study 1). Variables accepted in this model are hot ratio (0.0018181), Orel (- 
0.0000010), Dolol (-0.0000027), Ct(0.5 188709), 0re2(-0.0000069) and Constant(0.0249). 
Statistics of regression are presented in Table 5.1, model 5.4 (reg, method 2, Case Study 
1). Fig. 5.15 shows the graph of actual vs predicted carbon. It can be observed that 
R increased from 0.62 in model 5.4 (method2. Case Study 1) to 0.69 in model 5.4 (reg, 
method 2, Case Study 1) i.e., an improvement of 1 1.4%. The value of a also came down 
to 0.005 in model 5.4 (reg, method 2, Case Study 1) from 0.010 in model 5.4 (method 2, 
Case Study 1) i.e., a reduction of 50%. 

Model 5.4 (reg, method 2, Case Study 2) 

Model 5.4 (reg, method 2, Case Study 2) is regression model for model 5.4 (method 2, 
Case Study 2). Variables accepted in this model are Orel (-0.0000008), Dolol (-0.0000028), 
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Figure 5.13: Data from Plant 1, Case Study 1(273 heats): Actual vs Predicted for model 
5.4 (reg, method 1, Case Study 1), best fit lie is given by y=0.4469x+0.0322. Statistics are 
n=273, R=0.67, std. error=0.005. 
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Figure 5.14; Data from Plant 1, Case Study 2(156 heats): Actual vs Predicted for Model 
5.4 (reg, method 1, Case Study 2), best fit line is given by y=0.4333x+0.0331. Statistics 
are n=156, R=0.66, std. error=0.004. 
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Figure 5.15: Data from Plant 1, Case Study 1(273 heats): Actual vs Predicted for model 
5.4 (reg, method 2, Case Study 1), best fit line is given by y=0.0308x+0.4751. Statistics 
are n=273, R=0.69, std. error=0.005. 






ct(0.4660043) and Constant(0.0352). Statistics of regression are presented in Table 5.1 
(Model 5.4 (reg, method 2, Case Study 2)). Fig. 5.16 shows the graph of actual vs pre- 
dicted carbon. It can be observed from statistics of regression that R increased from 0.67 
in model 5.4 (method 2, Case Study 2) to 0.71 in model 5.4 (reg, method 2, Case Study 
2) i.e., an improvement of 5.5%. The value of a also came down to 0.004% in model 5.4 
(reg, method 2, Case Study 2) from 0.008% in model 5.4 (method 2, Case Study 2)i.e., a 
reduction of 50%. 

Model 5.5 (reg) 

Model 5.5 (reg) is regression model for model 5.5. Variables accepted in this model are 
Dolol (-0.0000043), To(0.0002185), Cf(0.5238060), Dolo2(-0.0000053), hlans2(0.0001689) 
and Constant(-0.3487). Statistics of regression are presented in Table 5.1, model 5.5 (reg). 
Fig. 5.17 shows the graph of actual vs predicted carbon. It can be observed from statistics 
of regression that R value increased from 0.37 in model 5.5 to 0.74 in model 5.5 (reg) i.e., 
an increase of 101%. The value of a also came down to 0.004 in model 5.5 (reg) from 
0.009 in model 5.5 i.e., a reduction of 51.7%. For a typical heat with dolol = 413, T(,= 
1607, dolo2 = 1501, hlans2 = 187, Ct = 0.068 

4 = -0.3487 + 0.5238060(0.068) - 0.0000043(413) + 0.0002185(1607) - 
0.0000053(1501) + 0.0001689(187) 

= 0.060. 

The actual carbon for this heat is 0.050. 

5.2.1.7 Comparison of various models of carbon prediction for Plant 1 

On comparing models 5.1 and 5.2, it can be observed that R value increased from 0.59 in 
model 5.1 to 0.66 in model 5.2. But the value of a also increased from 0.007 in model 5.1 
to 0.025 model 5.2, by 73%. On comparing model 5.1 (reg) and model 5.2 (reg), latter is 
able to predict better in terms of R and a (see table 5.1). This may be due to the fact that 
Qgff has been neglected in model 5.1. Also, by incorporating the operational variables 
like Orel, dolol in model 5.2 is able to give a better prediction. 

On Comparing models 5.2 and 5.3 it would be observed that, prediction is poorer 
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Figure 5.16: Data from Plant 1, Case Study 2(156 heats): Actual vs Predicted carbon 
for model 5.4 (reg, method 2, Case Study 2), best fit line is given by y=0.4962x+0.0291. 
Statistics are n=156, R=0.71, std. error=0.004. 
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Figure 5.17: Data from Plant 1, Case Study 2(156 heats): Actual vs predicted for model 
5.5 (reg), best fit line is y=0.55x+0.0261. Statistics are n=156, R=0.74, std. error=0.004. 
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for model 5.3. By incorporating operational variables like orel, dolol, hlans2, 022 and 
sublance carbon c© in model 5.3, prediction improved (see Table 5 . 1 , model 5.2 (reg) and 
model 5.3 (reg)). 

It can be observed from results of GA-decarb models that predictions for Case Study 
2 are better than Case Study 1 . This is because data set of Case Study 1 contains some 
heats in which ore 2 was added instead of dolo 2 . On Comparing model 5.4 (method 1 , 
Case Study 2) and model 5.4 (method 2 , Case Study 2) it can be seen that latter is able 
to predict better. This is because rj is allowed to vary as a function of gas evolution in 
method 2. There is 8.34% improvement in R and 12.34% improvement in from model 
5.4 (method 1 , Case Study 1 ) to model 5.4 (methodl. Case Study 2 ). Among GA-decarb 
models (see Table 5.1), model 5.4 (reg, method 2 , Case Study 2 ) is giving best predictions 
with R=0.71 and a =0.004%. Carefully observing the operational variables used in model 
5.4 (reg, method 2 , Case Study 2 ) and 5.3 (reg), lance height is used in model 5.3 (reg), 
which is missing in model 5.4(reg, method 2, Case Study 2 ). The possible reason for lower 
R value in model 5.4 (reg, method 2, Case Study 2) when compared to that of model 5.3 
(reg) could be exclusion of lance height. The is also supported by the fact that stirring 
energy transfered to the metal varies with lance height\ and in GA-decarb models 
Qeff is assumed to be constant. This is also the case with the other GA-decarb models. 
Variation of rj and mass transfer coefficient (K') are discussed latter. 

It can be observed that 77 value obtained from model 5.4, Case Study 2 is 2.477 where 
as that from model 5.5 is 2.91. It conforms that 77 value does change during the process. It 
is further observed that prediction is better with model 5.5 (reg) than with model 5.4 (reg, 
method 2, Case Study 2), with R 0.74 and nearly same a value. Hence, it confirms that K' 
is proportional to stirring energy, which decreases as a function of time. 

5.2.2 Temperature prediction models 

Temperature models are given by eq. (3.15) 


T = + (5.8) 

Z=1 

^See appendix 1 
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where T is predicted temperature and p. ’s are coefficients of operaUonal variables AT 
These models are developed using multiple linear regression (MLR). Temperature predic- 
tion models are tuned for Case Study 2 (156 data sets). ^ 


Model 5.6 

Fig 5.18 shows graph of actual versus predicted temperature. From the statistics of regres- 
sion (Table 5.2) it can be observed that c = 5 “C and R=0.90 (0.9027).Thus, this model 
is able to correlate 81% of variation in the data. Operational variables used and their co- 
efficients are hot ratio( 1.6089489), r<,(0.8260166), 022(0.0349623), slag2(-0.0051884) 
dolo2(-0.0105275) and constant(265.4443). For a typical heat with hot ratio = 5 2098 
022 = 2464, To= 1607, slg2 = 1508, dolo2 = 1501 

T = 265.4443 -i- 1.6089489(5.2098) + 0.8260166(1607) + 0.349623(2464) - 
0.0051884(1508) - 0.0105275(1501) 

= 16630c. 

The actual temperature for this heat is 1662 ^C. 


5.2.3 Dissolved oxygen prediction models 

Dissolved oxygen models are given by eq. (3.15) 


n— 1 

[O] = Mo + ^ Ml • (5 9) 

t=l 

where ’s are coefficients of operational variables X,. These models are developed using 
MLR. Dissolved oxygen models are tuned for Case Study 2 (156 heats). 

Model 5.7 

Table 5.3, model 5.7 shows the statistics of regression for actual versus predicted dissolved 
oxygen. Fig 5.19 shows graph of actual versus predicted dissolved oxygen. From statis- 
tics of regression it can be observed that cr=37 ppm and R=0.83 (0.8267). Thus this model 
explains 68.3% variation in data. Operational variables used are basicity(12.3357568), 


86 



CO 

O 3 

CO 


S 

i 


Kgs - N 

H S « § 'o 

bb &-C ^ 

<D ST cd ^ 

T 3 > 


^ "S 

I ° N « 

e 3 O .§ 


JL (2 *13 ^ 

3 4) 25 S o iS 

0- 2 ® ^ Q Tl w 


II , *-H 

e>h c.a 

w 


87 



Predicted Temperaure 


1660 


1650 


>oO^.C 8> 


o o 

/sOa.O 




1640 


.•• 0 £> 


O O 

1630 ^ 

1630 1640 


1650 1660 1670 1680 1690 17 ( 

Actual Temperature 


Figure 5.18: Data from Plant 1, Case Study 2(156 heats): Actual vs Predicted tempera- 
ture for model 5.6, best fit line is y=0.8146x+307.76. Statistics are n=156, R=0.90, std. 
eiTor=5. 



Co("895.1580774), 022(0.2812597) and constant( 123. 6482). For a typical heat with ba- 
sicity=5.443, =0.326, 022=2464 

[O] = 123.6482+12.3357568(5.443)-895.1580774(0.326)+0.2812597(2464) 

=592 ppm. 

The actual dissolved oxygen content for this heat is 413 ppm. 

In order to improve the prediction, lance height during the second blow (hlans2) and 
end point temperature (T) are also considered as independent variables. The results are as 
presented in model 5.8. 

Model 5.8 

From statistics of regression (Table 5.3, model 5.8) it can be observed that a=34 ppm 
and R=0.87 (0.8724). Fig 5.20 shows graph of actual versus predicted dissolved oxy- 
gen. Thus this model explains 76% variation in data. It can be observed from model 
5.7 and model 5.8 that prediction improved in terms of both R and a. lliere is an im- 
provement of 5.5%(approx) in R and 7.6% in o in model 5.8 from model 5.7. Opera- 
tional variables used and their coefficients are basicity(12.9760397), c<,(-875.4523083), 
022(0.1364417), 7;(-2.9485933), dolo2(0.0428712), hlans2(-0.4622155), T(3.3075463) 
and constant(-300.6208). For a typical heat with basicity = 5.443, To=- 1607, 022 = 2464, 
ci= 0.517, dolo2 = 1501, hlans2 = 187, T = 1663 (see model 5.6) 

[O] = -300.6208 + 12.9760397(5.443) - 2.9485933(1607) - 875.45230'83(0.517) + 
0.1364417(2464) + 0.0428712(1501) - 0.4622155(187) + 3.3075463(1663) 

= 394 ppm. 

The actual dissolved oxygen content for this heat is 413 ppm. 

5.2.4 Phosphorus prediction models 

Phosphorus models have been discussed in section 2.3.5. Detailed study was done first on 
Plant 1 and the best model is then tuned for Plant 3 data. In case of Plant 2, phosphorus 
content of heats was not available. 
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Figure 5.19: Data from Plant 1, Case Study 2(156 heats): Actual vs predicted [O] for mod- 
el 5.7, best fit line is y=0.6786x+148.1316. Statistics are n=156, R=0.83, std. error=37. 
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Figure 5.20: Data from Plant 1, Case Study 2(156 heats): Actual vs predicted [O] for 
model 5.8, best fit line is y=0.761x+101.255. Statistics are n=156, R=0.87, std, error=34. 




5.2.4.1 Kinetic model of mass transfer 


Kinetic model of mass transfer is given by eq. (3.21) 

“ (5.10) 

where ’s are coefficients of operational variables X, . These models are developed using 
MLR. 

Model 5.9 

In this model, equilibrium constant of phosphorus is taken as 2. From statistics of re- 
gression (Table 5.4, model 5.9) it can be observed that (t= 0.004 (0.00352)% and R=0.60 
(0.6026). This model is able to explain 36.3% of the variation. Operational variables 
used and their coefficients are basicity(-0.0004353), dolo2(-0.0000018), Co(-0.0099315), 
hlans2(0.0000365), 022(0.0000030) and constant(-0.0126). Fig. 5.21 shows graph of ac- 
tual versus predicted phosphorus. For a heat with basicity=5.443, hlans2=187, 022=2464, 
Co=0.326, dolo2=1501we obtain by substituting in eq. (5.10) 

- ln( )=-0.0 1 26-0.0004353(5. 443)-0.00000 1 8(1 501 )-0.00993 1 5(0.326)+ 

0.0000365(187)+0.0000030(2464)=-6.69xl0-® 

^ = 1.0067 
Pt =2+(0.027-2)( 1.0067) 

=0.014 

Actual end point phosphorus value for this particular heat is 0.008. 

Model 5.10 

From statistics of regression on actual versus predicted phosphoms (Table 5.4, model 5.10) 
it can be observed that <7=0.003% and R=0.87. Fig 5.22 shows graph of actual versus pre- 
dicted phosphorus. This model is able to explain 75.7% of the variation. It is observed 
from model 5.9 and 5.10 that by excluding the equilibrium factor of phosphoms from 
the dependent variable, prediction improved from 0.60 in model 5.9 to 0.87 in model 
5.10 in terms of R and cr decreased to 0.0034 in model 5.10 from 0.0035 in model 5.9. 
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Table 5.4 Summary of results obtained for phosphorus prediction for Plant 1 
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Predicted Phosphorus (mass%) 
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ngure 5.21: Data from Plant 1, Case Study 2(156 heats): Act 
•us for model 5.9, best fit line is y=1.3253x+0.0011. Statistic 
5rror=0.003. Each square may represent more than one heat. 




Operational variables used and their coefficients are svol(0.0000181), Orel (-0.0000140), 
slg2(0.0000968), hlans2(-0.0013645), dolo2(0.0001475), Mn2(-1.8107433) and constan- 
t(1.1448). For a typical heat with svol= 17208.47, orel=3994, slg2=1508, hlans2=187, 
dolo2=1501, Mn2=0.174 we obtain 

“ •l'^^+0-000018(17208.47)-0.0000140(3994)+0.0000968(1508)- 

0.0013645(187)+0.0001475(1501)-1.8107433(0.174)=1.1975 
pt =0.027(0.302) = 0.008. 

Actual phosphorus for this heat is 0.008. 

5.2.4.2 Regression model 
Regression model is given by eq. (3.22) 


n-1 

p[ = fXo + l^iPo + /Wi • (5.11) 

1=2 

where, pj is the end point phosphorus value predicted from MLR, po is sublance phospho- 
rus and Pi ’s are coefficients (to be determined by MLR) of operational variables Xi 

Model 5.11 

From statistics of regression on actual versus predicted phosphoms (Table 5.4, model 5.1 1 
)it can be observed that (7=0.0007% and R=0.93. Fig 5.23 shows graph of actual ver- 
sus predicted phosphorus. This model is able to explain 85.6% of the variation. Op- 
erational variables used and their coefficients are svol(-0.0000002), htr * po(0.0032876), 
slg2(0.0000968), hlans2(-0.0013645),po(0.1063842), Mn2(0.0445678), T(0.0000339) and 
constant(1.1448). For a typical heat with slag mass = 17208.47, T = 1663 (see model 5.6), 
htr*p = 0.2917, slg2 = 1508, Mn2 =.172 , hlans2 = 187, Po = 0.027 

p't = -0.0533 - 0.0000002(17208.47) + 0.0032876(0.2917) + 0.1063842(0.027) - 
0.0000004(1508) - 0.0000081(187) + 0.0000339(1663) + 0.0445678(0.172) 

= 0.009. 

Actual phosphorus for this heat is 0.008. 
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5.2.4.3 Discussion on phosphorus prediction results 

It is observed that the best model for phosphorus prediction is model 5.11. Hence, this 
model is applied to Plant 3. In case of Plant 2, we did not have the data for measurement 
of phosphorus. Hence for only Plant 3 phosphorus model is developed. 

5.2.5 Manganese prediction models 

Manganese models developed in this work are discussed in section 3.9.2. 

5.2.5.1 Kinetic model of mass transfer 
Kinetic model of mass transfer is given by eq. (3.26) 


Mnj — Mne 
Mrio — Mue 


n~l 

f^o + ^ fJn. ■ 
1=1 


where fj,^ ’s are coefficients of operational variables X^. 


(5.12) 


Model 5.12 

In this model, equilibrium constant of manganese Mrie is taken as 0.005. Table 5.5, mod- 
el 5.12 shows the statistics of regression for actual versus predicted manganese. Fig 5.24 
shows graph of actual versus predicted manganese. From statistics of regression it can 
be observed that cr=0.007 (0.00699)% and R=0.96 (0.9608). This model is able to ex- 
plain 92.2% of the variation. Operational variables used and their coefficients are To(- 
0.0027084), Co(-0.8401064), 022 (0.0002366), hlans2 (-0.0011956), dolo2 (0.0001237), 
T (0.0035221), pt(-43.9292138) and constant(-0.3590). For a typical heat with To=1607, 
Co=0.326, 022=2464, hlans2=187, dolo2=1501, T=1662,pt=0.008 we obtain 
_ In MaL;:»=.o.359.0.0027084(1607)-0.8401064(0.326)+0.0002366(2464)- 

TIq — 'U uUo 

0.001 1956(187)4- 0.0001237(1501)+0.003522(1662)-43.9292138(0.008)=1.062 
Mnt =(0.517-0.005)(0.3457)-K).005 
= 0.182%. 

Actual manganese for this heat is 0.174 %. 
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5.2.4.3 Discussion on phosphorus prediction results 

It is observed that the best model for phosphorus prediction is model 5.11. Hence, this 
model is applied to Plant 3. In case of Plant 2, we did not have the data for measurement 
of phosphorus. Hence for only Plant 3 phosphorus model is developed. 


5.2.5 Manganese prediction models 

Manganese models developed in this work are discussed in section 3.9.2. 


5.2.5.1 Kinetic model of mass transfer 
Kinetic model of mass transfer is given by eq. (3.26) 


Mrit — Mrig 
Muo - Mn, 


~ d" ^ ] Mt ■ 
1=1 


where fXt ’s are coefficients of operational variables X,. 


(5.12) 


Model 5.12 

In this model, equilibrium constant of manganese Mrig is taken as 0.005. Table 5.5, mod- 
el 5.12 shows the statistics of regression for actual versus predicted manganese. Fig 5.24 
shows graph of actual versus predicted manganese. From statistics of regression it can 
be observed that cr=0.007 (0.00699)% and R=0.96 (0.9608). This model is able to ex- 
plain 92.2% of the variation. Operational variables used and their coefficients are Td- 
0.0027084), Co(-0.8401064), 022 (0.0002366), hlans2 (-0.0011956), dolo2 (0.0001237), 
T (0.0035221), pt (-43. 9292 138) and constant(-0.3590). For a typical heat with To=1607, 
Co=0.326, 022=2464, hlans2=187, dolo2=1501, T=1662, pt=0.008 we obtain 
- In f2t=|^=.0.359-0.0027084(1607)-0.8401064(0.326)+0.0002366(2464)- 
0.001 1956(187)+ 0.0001237(1501)+0.003522(1662)-43.9292138(0.008)=1.062 
Mnt =(0.517-0.005)(0.3457)+0.005 
= 0.182%. 

Actual manganese for this heat is 0.174 %. 


99 



a TJ 

U 3 

03 


03 

G 3 nd 

u 3 

03 


03 ^ 

03 TJ 

u 3 


On ^ 

S 00 ^ 
o On 

^ q d 

9 o II 

a !1“ 


§5“SS??S! 

b ir M jl ^ - n 


»n ^ JT; 
^ On o r- 
O o ^ 

O II O 00 

jl C^ 1} 1\ 


Sm 4-J 

O G 

-S ii - -- 
S G -J S 

feb *c kS 
P Sr c3 ^ d 

^ 'O >• ^ d' 

4 -( c 

o 


•a a 
.2 S i> 

s I •§ fi* 00 

&l^i 

Q 


+-J I 

S S- o 

Son 

5 2 *s 

c3 53 .S 

D. 


■S h! 


c 5 | o u! 


r-. ^ 

T3 o 

T 3 *+3 
O cd 

<D 


i 

II 

" ill 


: 3 . * 
+ ^ 


c<» 

+ = 3 . 


100 



Predicted Manganese (mass%) 





5.2.5.2 Regression model 
Regression model is given by eq. (3.27) 


Mn[ = fJ.o + iJ-iMrio + Ea**-^* (5-13) 

1=2 

where, Mn* is the end point manganese value predicted from MLR, Muo is the sublance 
manganese and ^x ^ ’s are coefficients of operational variables Xj 

Model 5.13 

From statistics of regression on actual versus predicted manganese (Table 5.5, model 5.13 
) it can be observed that a=0.001% and R=0.94. Fig 5.25 shows graph of actual versus 
predicted manganese. This model is able to explain 88.4% of the variation. Operational 
variables used and their coefficients are (0.00030 10), Co(0. 1810235), 022(-0.0000522), 
Orel (-0.0000012), dolo2(-0.000014), Mrio (0.4679948) and constant(-0.4030). For a typ- 
ical heat with orel = 3994, To= 1607, c,, = 0.517, Muo = 0.326, 022 = 2464, dolo2 = 
1501 

Writ = -0.4030 + 0.4679948 (0.0326) - 0.0000012(3994) + 0.0003010(1607) + 
0.1810235(0.517) - 0.0000522(2464) - 0.0000140(1501) 

= 0.172. 

The actual manganese for this heat is 0.174. 

Model 5.14 

From statistics of regression on actual versus predicted manganese (Table 5.5, model 5.14 
) it can be observed that cr=0.0065% and R=0.96. Fig 5.26 shows graph of actual versus 
predicted manganese. This model is able to explain 92.2% of the variation. Operational 
variables used are pt(4.993776), Co(0. 10767 12), 022(-0.0000408), hlans2(0.0001395), 
dolo2(-0.0000106), Mrio (0.4075651) and constant(0.0233). For a typical heat with c<, 
=0.326, 022=2464, hlans2=187, dolo2=1501, Mno=0.517, pt=0.008 
Mn\ = 0.0233+0.1076712(0.326)-0.0000408(24640+0.0001395(187)- 
0.0000106(1501)+0.4075651(0.5170+4.993776(0.008) 
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=0.218 

Actual manganese for this heat is 0.174. 

5.2.5.3 Discussion on manganese prediction results 

It is observed that prediction from model 5.14 is better than model 5.13 and model 5.12. 
It can be observed from results of model 5.12 and model 5.14 that although the R value is 
same^ cT value is less in latter model (0.007, 0.0065). Hence this model is further tuned on 
Plant 3 data. 

5.3 Results of Plant 2 

5.3.1 End point carbon prediction models 

5.3.1.1 GA-decarb models (model 5.4) 

GA-decarb models tuned in this data in this section. 

Model 5.4 (method 1) 

Eq. (5.4) was optimized for 340 data sets. Optimized values of ai and rj corresponding 
to minimum SSE are 0.0084 and 0.2910. Table 5.6, model 5.4 shows the statistics of 
regression performed on actual vs predicted carbon for Plant 2. From the statistics of 
regression it can be observed that a=0.012 (0.0117)% and R= 0.61 (0.6129). Fig. 5.27 
shows actual vs predicted graph. 

Model 5.4 (method 2) 

Eq. (5.5) was optimized for 340 data sets. Optimized values of ai and Oa corresponding to 
minimum SSE are 0.017 and 0.0014. Table 5.6 (method 5.4, method 2) shows the statistics 
of regression performed on actual vs predicted carbon. It can be observed from the table 
that <7=0.014% and R=0.71 (0.7060). Fig. 5.28 shows graph of actual vs predicted carbon. 
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Figure 5.26: Data from Plant 1, Case Study 2(156 heats): Actual vs predicted manganese 
for model 5.14, best fit line is y=0.9237x+0.0115. Statistics are n=156, R=0.96, std. er- 
ror=0.0065. 
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Figure 5.27: Data from Plant 2 (340 heats): Actual vs predicted carbon for model 5.4 
(method 1), best fit line is y=0.4739x+0.0249. Statistics are n=340, R=0.61, std. er- 

ror=0.012. 
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Table 5.6 Summary of results obtained for carbon prediction for Plant 2 
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Figure 5.28: Data from Plant 2 (340 heats): Actual vs predicted carbon for model 5.4 
(method 2), best fit line is y=0.7433x+0.0145. Statistics are n=340, R=0.71, std. er- 
ror=0.014. 
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5.3.1.2 Regression models 
Model 5.4 (reg, method 1) 

Model 5.4 (reg, method 1) is regression model for model 5.4 (method 1). Variables ac- 
cepted in this model are Scrap (-0.0000006), Ct(0.6997620), chill scrap(0.0000029) and 
Constant(0.0311). Statistics of regression are presented in Table 5.6 (method 5.4 (reg, 
method 1)). Fig. 5.29 shows the graph of actual vs predicted carbon. It can be observed 
that R increased from 0.61 in model 5.4 (method l)to 0.64 in model 5.4 (reg, method 1) 
i.e., an improvement of 4.2%. a also came down to 0.010 in model 5.4 (reg, method 1) 
from 0.012 in model 5.4 ( method 1) i.e., a reduction of 13%. 

Model 5.4 (reg, method 2) 

Model 5.4 (reg, method 2) is regression model for model 5.4 (method 2). Variables accept- 
ed in this model are Scrap (-0.0000004), q( 0.5904476), chill scrap(0.0000017) and Con- 
stant(0.0282). Statistics of regression are presented in Table 5.6 (model 5.4 (reg, method 
2)). Fig. 5.30 shows the graph of actual vs predicted carbon. It can be observed that R 
value decreased slightly from 0.7060 in model 5.4 ( method 2) to 0.7050 in model 5.4 (reg, 
method 2). a also came down to 0.010 in model 5.4 (reg, method 2) from 0.014 in model 
5.4 ( method 2) i.e., a reduction of 30%. 

Model 5.3 (reg) 

Model 5.3 (reg) is developed based on 321 heats (with end point carbon (c*) in the range 
0.02-0.075) of Plant 2. First linear regression was performed to fit best fit line for model 
5.3. Fig. 5.3 1 shows actual vs predicted carbon graph. But after observing the trend of re- 
sults in Fig. 5.31, polynomial of 2nd order was tried. Then again multiple linear reression 
was performed to obtain a model with statistics R=0.81 and a=0.01 1 . After observing the 
parameters accepted, MLR is directly performed on the variables. The model so evolved 
is the one discussed here. 

This model based on polynomial fit of 2nd order to the Linear model. Variables ac- 
cepted in this model are 022(0.0000000186), 022(-0.0000421756), ro(0.0000688766). 
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Figure 5.30; Data from Plant 2 (340 heats): Actual vs Predicted carbon for model 5.4 
(reg, method 2), best fit line is y=0.5064x+0.0201. Statisatics are n=340, R=0.71, std. 
error=0.010. 
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Figure 5.31; Data from Plant 2 (321 heats): Actual vs predicted carbon for model 
Best fit line is y=0.5077x+0.0215. Statistics are n=339, R=0.71, std. error=0.010. 



flux(-0.0000007942), c2(0.06767 18982), On • Co(-0.00008 18279), Co(0. 1238672564) and 
constant(-0.05). Statistics of regression are presented in Table 5.6 (method 5.3 (reg)). Fig. 
5.32 shows the graph of actual vs predicted carbon. It can be observed from the statistic- 
s that (7=0.006 (0.00554) and R=0.74 (0.7414). For a typical heat with 0222=4524129, 
c2=0.202, 022 * Co =957.150, 022=2127, Co=0.450, To=1534, flux=12450 

c;=-0.05+0.0000000186(4524129)-0.0000421756(2127)+0.0000688766(1534)- 

0.0000007942(12450)+0.0676718982(0.202)-0.0000818279(957.150)+ 

0.1238672564(0.450) 

=0.028. 

Actual end point carbon value for this heat is 0.029. 

5.3.1.3 Discussion on carbon prediction results 

On Comparing results of prediction by model 5.4 (method 1) and model 5.4 (method 2), it 
can be observed that the latter model is giving better results with R of 0.7 1 and a of 0.014. 
This is 15% improvement in R over that of model 5.4 (method 1). This improvement is 
because rj is allowed to vary in method 2 of GA-decarb models. After performing MLR 
on the results obtained from GA-decarb models, value of cr decreased to 0.010 in model 
5.4 (reg, method2) from 0.014 in model 5.4(method 2). 

The best prediction for Plant 2 is given by model 5.3 (reg) with R of 0.7414 and a of 
0.00554%. This cr is 43% lower than the one obtained from model 5.4 (reg, method 2). 

5.3.2 Temperature prediction results 

Model 5.6 

From the statistics of regression on actual versus predicted temperature(Table 5.7) it can 
be observed that a = 8° C and R=0.72 (0.7241). Fig. 5.33 shows graph of actual versus 
predicted temperature. This model is able to correlate 52% of variation in the data. Opera- 
tional variables used are 022(0.0230991), ro(0.58), Co( 14. 127787), flux(-0.0009627) and 
constant(7 14.4507). For a typical heat with 022=2127, Cq =0.450, To=1534, flux=12450 
T=714.4507+0.0230991(2127)+0.58(1534)-hl4.127787(0.45)-0.0009627(12450) 
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Figure 5.32: Data from Plant 2 (321 heats): Actual vs predicted carbon for model 5.3 
(reg). Best fit line is y=0.5498x+0.0185. Statistics are n=321, R=0.74, std. errop=0.005. 



=1648. 

Actual temperature for this heat is 1670. 

5.3.3 Dissolved oxygen prediction results 

Model 5.8 

From the statistics of regression on actual versus predicted dissolved oxygen (Table 5.8 
) it can be observed that ct= 139 ppm and R=0.82 (0.8181). Fig. 5.34 shows graph of 
actual versus predicted dissolved oxygen. This model is able to explain 67% of the varia- 
tion. Operational variables used are T(8.7309033), Co(-1542.8043187), 022(0.2338949), 
To(-7. 1812861) and constant(-2002.4246). For a typical heat with T=1670, c^, =0.450, 
022=2127, To =1534 

[0]=-2002.4246+8.7309033(1670)-1542.8043187(0.45)+0.2338949(2127)- 

7.1812861(1534) 

=1365 ppm. 

Actual dissolved oxygen for this heat is 1432 ppm. 


5.4 Results of Plant 3 

The models developed for this Plant are tuned for a total of 297 heats . 

5.4.1 End point carbon prediction models 

5.4.1.1 GA-decarb models (model 5.4) 

Ga-decarb models are presented in section 5.2. 1.4. The result for method 1 of GA-decarb 
models is very poor and hence not presented here. 


Model 5.4 (method 2) 

This is GA-decab model developed as explained in eq. (5.5). Model is optimized for 297 
data sets. Optimized values of oji and corresponding to minimum SSE are 0.0195 and 
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Figure 5.33: Data from Plant 2 (321 heats): Actual vs Predicted temperature for Model 
5.6, best fit line is y=0.5269x+786.2107. Statistics are n=321, R=0.72, std. errors:8.0. 
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Figure 5.34: Data from Plant 2 (321 heats): Actual vs predicted [0] for model 5.8, best fit 
line is y=0.6695x+296.0932. Statistics are n=321, R=0.82, std. error=140. 
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0.00175. Table 5.9 shows the statistics of regression performed on actual vs predicted 
carbon for paint 3. From the statistics of regression it can be observed that cr=0.016 and 
R=0.21 (see Table 5.9, model 5.4 (method 2)). 

5.4.1.2 Regression models 
Model 5.4 (reg, method 2) 

Model 5.4 (reg, method 2) is regression model for model 5.4 (method 2). Variables 
accepted in this model are basicity (-O.OOIXXIIO), Ct(0.3 133347), To(-0.0001404), 022 
(0.0000152) and constant (0.2630). Statistics of regression are presented in Table 5.9 
(model 5.4 (reg, method 2)). Fig. 5.35 shows the graph of actual vs predicted carbon. It 
can be observed from statistics of regression that R value increased from 0.21 in model 5.4 
(method 2)to 0.44 in model 5.4 (reg, method 2) i.e., an increase of 107%. o also reduced 
to 0.004 in model 5.4 (reg, method 2) from 0.016 in model 5.4 (method 2) of Plant 3 i.e., 
a reduction of 76.3%. 

Model 5.3 (reg) 

Model 5.3 (reg) is tuned for 297 heats of Plant 3. Variables accepted in this model 
are basicity(-O.OOOOll), ro(-0.0001478), Co(0.0232972), 022(-0.0000227) and constan- 
t(0.3038). Statistics of regression performed on actual and predicted carbon are presented 
in Table 5.9 (model 5.3 (reg)). From statistics of regression it can be observed that <7=0.003 
(0.0034)% and R=0.39 (0.3783). Fig. 5.36 shows the graph of actual vs predicted carbon. 
For a typical heat with basicity=3.85, ro=1629, Co =0.549, 022=740 

Ct = 0.3038-0.00001 1(3.85)-0.0001478(1629)+0.0232972(0.549)-0.0000227(740) 
=0.059. 

Actual end point carbon for this heat is 0.050. 

5.4.1.3 Discussion on carbon prediction results 

From prediction statistics of model 5.4 (reg, method 2) and model 5.3 (reg) it can be 
observed that R value decreased from 0.44 in model 5.4 (reg, method 2) to 0.38 in model 
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Figure 5.35: Data from Plant 3 (297 heats): Actual vs predicted carbon for model 5.4 (reg, 
method2), best fit line is y=0.1929x+0.041. Statistics are n=297, R=0.44, std. error=0.(X)4. 
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Figure 5.36: Data from Plant 3 (297 heats): Actual vs predicted carbon for model 5.3(reg), 
best fit line is y=0.1435x+0.0441. Statistics are n=297, R=0.39, std. error^O.003. 
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5.3 (reg). But a value improved from 0.004 in model 5.4 (reg, method 2) to 0.003 in model 
5.3 (reg). 


5.4.2 Temperature prediction results 

Model 5.6 

From statistics of regression on actual versus predicted temperature (Table 5.10) it can be 
observed that a = 9 C and R=0.73 (0.7310). Fig. 5.37 shows graph of actual versus 
predicted temperature. This model is able to correlate 59.6% of variation in the data. 
Operational variables used are raw dolo2(-0.0289692), To(0.8 109563), 022(0.0982545), 
hlans2(0.6512143) and constant(215.8787). For a typical heat with raw dolo2=650, To 
=1629, 022=1AQ, hlans2=148 

T = 215.8787-0.0289692(650)+0.8109563(1629)+0.0982545(740)+ 

0.6512143(148) 

=1687. 

Actual temperature for this heat is 1702. 

5.4.3 Dissolved oxygen prediction results 

Oxygen models are developed for 297 heats. 

Model 5.8 

From statistics of regression on actual versus predicted dissolved oxygen (Table 5.11 ) 
it can be observed that o-ll ppm and R=0.60 (0.5978). Fig. 5.38 shows graph of actual 
versus predicted dissolved oxygen. This model is able to explain 35% of the variation. Op- 
erational variables used are basicity(0.0175579), llife(0.2249089), raw dolo2(0.0847864), 
ro(-0.9097412), T(5.2328732) and constant(-6928.4438). For a typical heat with basici- 
ty=3.85, llife=6, raw dolo2=650. To =1629, T=1702 

{0] = -6928.4438+0.0175579(3.85)+0.2249089(6)+0.0847864(650)-0.9097412(1629)+ 
5.2328732(1702) 

=552. 
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Figure 5.37: Data from Plant 3 (297 heats); Actual vs Predicted temperature for model 
5.6, best fit line is y=0.5351x+785.333. Statistics are n=297, R=0.73, std. error=9. 



Actual dissolved oxygen for this heat is 680. 


5.4.4 Phosphorus prediction results 

Phosphorus models are developed for 297 heats. 

Model 5.11 

From statistics of regression on actual versus predicted phosphorus (Table 5.12) it can be 
observed that cr=0.001% and R=0.81 (0.8113). Fig. 5.39 shows graph of actual versus 
predicted phosphorus. Using this model 75% of variation can be explained. Operational 
variables used are svol(-0.0000007). To (-0.0000237), Co(-0.00 16602), 022(-0.0000027), 

PI (0.2280 134), Mn2(0.0237083), T(0.0000471) and constant (-0.0357 ). For a typical 
heat with svol=2871 .444, 022=740, To =1629, Co =0.549, Pl=0.018, Mn2=0.120, T=1702 
p' = -0.0357-0.0000007(2871 .444)-0.0000237(1629)-0.0016602(0.549)-0.0000027(740)+ 
0.2280 1 34(0.01 8)+0.0237083(0. 120)+0.000047 1 ( 1 702) 

=0.008 

Actual phosphorus for this heat is 0.009. 

5.4.5 Manganese prediction results 

5.4.5.1 Regression model 

Regression model is given by eq. (3.27) 

n—l 

Mn't = iio + + Y^^r-X, (5.14) 

2=2 

where, Mn[ is the end point manganese value predicted from MLR, Muo is the sublance 
manganese and fii ’s are coefficients of operational variables Xi 


Model 5.13 

From statistics of regression on actual versus predicted manganese (Table 5.13) it can be 
observed that a=0.012% and R=0.83. Fig. 5.40 shows graph of actual versus predicted 
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Figure 5.38: Data from Plant 3 (297 heats); Actual vs predicted [O] for model 5.8, best fit 
line is y=0.3571x+387.1967. Statistics are n=297, R=0.60, std. error=77. 
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manganese. Variables accepted in this model are Mno(0.4880717), 022(-0.0000294), 
raw dolo2(-0.0000091), hlans2(0.0006755) and constant(-0.0436). For a typical heat with 
Muo =0.19, 022=740, raw dolo2=650, hlans2=148 

M nj =-0.0436+0.48807 1 7(0. 19)-0.0000294(740)-0.000009 1 (650)+0.0006755( 148) 
= 0 . 121 . 

Actual manganese for this heat is 0.120. 

5.5 Comparision of models developed for Plant 1, Plant 2 
and Plant 3 

In this section the best models tuned for Plants 1, 2 and 3 are compared to estimate the 
performance of each plant with respect to another. 

5.5.1 Discussion on carbon prediction models 

Models which are in common to all three plants are GA-decarb models (model 5.4) and 
model 5.3(reg). From model 5.4( reg, methodl. Case Study 2 ) of Plant 1 and model 5.4 
(reg, method 1) of Plant 2, it can be observed that the former is giving better prediction 
with cr=0.004%C. Among model 5.4 (reg, method 2, Case Study 2) of Plant 1, model 5.4 
(reg, method 2) of Plant 2 and model 5.4 (reg, method 2) of Plant 3, Plant 1 has lowest 
value of a (0.004), giving the best prediction. 

From model 5.3 (reg) of Plants 1, 2 and 3 it can be observed that standard error a is 
lowest (0.0034) in case of Plant 3. Best prediction is given by model 5.3 (reg) of Plant 
1 with R value of 0.7617 and cr of 0.00418. Observing the parameters used in model 5.3 
(reg), lance height is selected in Plant 1 and not in Plant 3. In case of Plant 2, lance height 
is not considered as it is kept constant for all the heats. Summary of best results for carbon 
prediction are summarized in Table 5.14, appendix 2, appendix 3 and appendix 4. 
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Figure 5.40: Data from Plant 3 (297 heats): Actual vs predicted manganese for model 
5.13, best fit line is y=0.6848x+0.0456. Statistics are n=297, R=0.83, std. error=0.012. 
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F=70.494 



5.5.2 Discussion of temperature prediction models 

Frorn temperature models developed for Plants 1, 2 and 3 it can be observed that model 
5.6 of Plant 1 is giving best prediction with R of 0.9027 and a of 5.22 ^C. The error of 
prediction in case of Plants 2 and 3 is 8 and 9 , respectively. Summary of best 

results for temperature prediction are summarized in Table 5.15, appendix 2, appendix 3 
and appendix 4. 

5.5.3 Discussion of dissolved oxygen prediction models 

From dissolved oxygen models developed for Plants 1, 2 and 3 it can be observed that 
model 5.8 of Plant 1 is giving best prediction with R of 0.87 and a of 34 ppm. Although R 
value is better for model 5.8 of Plant 2, a is better for Plant 3. The accuracy of prediction 
in case of Plants 2 and 3 are 139 ppm and 77 ppm , respectively. Summary of best results 
for dissolved oxygen prediction are presented in Table 5.16, appendix 2, appendix 3 and 
appendix 4. 

5.5.4 Discussion of phosphorus prediction models 

No data was available for phosphorus prediction in case of Plant 2. Hence the models 
developed for Plants 1 and 3 are compared in this section. Comparing model 5.11 of 
Plants 1 and 3, Plant 1 is giving better prediction with R of 0.93 and a of 0.0007%. In case 
of Plant 3, accuracy of phosphoras predicition is 0.001%. The best results for phosphoras 
prediction are summarized in Table 5.17, appendix 2, appendix 3 and appendix 4. 


5.5.5 Discussion of manganese prediction models 

No data was available for manganese prediction for the case of Plant 2. Hence the models 
developed for Plants 1 and 3 are compared in this section. Comparing model 5 . 1 3 of Plants 
1 and 3, Plant 1 is giving better prediction with R of 0.96 and cr of 0.0065%. In case of 
Plant 3, accuracy of manganese predicition is 0.012%. The best results for manganese 
prediction are summarized in Table 5.18, appendix 2, appendix 3 and appendix 4. 
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F=165.1967 
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5.6 Evaluation of capacity mass transfer coefficient {K') 

Capacity mass transfer coefficienthas been assumed , as suggested literature[13], to 
depend on gas flow rate as follows 

( dc 

--^Pi + Qeffj (5.15) 

where 77 is exponential factor. In the present work value of K has been determined by 
two methods by using GA-decarb models for Plant 1 , Plant 2 and Plant 3 . The results are 
discussed for each plant followed by comparison of all three plants. 

5.6.1 Capacity Mass transfer coefficient for Plant 1 

GA-decarb model (Method 1) 

In this method rj in eq. (5.15) is assumed to be constant through out the process. Fig. 5.41 
shows the variation of mass transfer coefficient with time. From model 5.4 (method 1, 
Case Study 1) and model 5.4 (method 1, Case Study 2) it can be observed that while oi 
values are very close to each other, the rj values are different in both cases (although they 
are from same Plant) as seen from Table 5.19. This is because the data in Case Study 1 
has some heats in which ore2 has also been added during second blow as coolant. Hence, 
there is combined effect of ore2 and dolo2 in optimizing the rj value in Case Study 1, but 
in Case Study 2 only dolo2 plays role. Because of this difference in rj, there is difference 
in the /(''values as observed from Fig- 5.41(a) and 5.41(b) and the values in Fig. 5.41 
(a) are slightly lower than in Fig. 5.41(b) in first 100 seconds. This may be attributed to 
higher slag volume in the former case due to ore addition. 

GA-decarb model (Method 2) 

In this method 77 is made to vary as a function of gas flow rate as 
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Capacity mass transfer coefficient (K , per sec ) 



(a) Co = 0.369 and Ct=0.044 (b) Co = 0.369 and Ct=0.044 
(c) Co = 0.365 and Ct=0.043 


^ (K') VS time for method 1 of GA-decarb 

Figure 5.41: Capacity mass transfer coetticieni ; vb uiuc lui 

models. 
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Variation of 77 as a function of time is shown in Fig. 5.42. This graph is drawn for heat 
with Co=0.369 and Cf=0.044. From Fig 5.42 (a) and 5.42 (b) it can be observed that the 
difference in 77 value between Case Study 1 and Case Study 2 is narrowing down towards 
the end of blow and the values in Case Study 1 are lower than in Case Study 2. Fig. 5.43 

(a) shows the variation of mass transfer coefficient, K' with time. From Fig. 5.43(a) and 

(b) it can be observed that K'is. almost the same for model 5.4 (method 2 , Case Study 
1 ) and model 5.4 (method 2 , Case Study 2 ). It can also be observed that mass transfer 
coefficient decreases almost linearly with time but at a very slow rate. This is because gas 
evolution rate decreases with time. 


5.6.2 Capacity Mass transfer coefficient for Plant 2 

Method 1 

Fig. 5.41 (c) shows the variation of mass transfer coefficient K' with time for model 5.4 
(method 1) of Plant 2. 

Method 2 

Variation of 77 with time is shown in Fig. 5.42 (b) and Fig. 5.43 (c) shows the variation 
of mass transfer coefficient {K') with time. It can be observed that, similar to plant 1, K' 
almost linearly decreases with time. 

5.6.3 Capacity Mass transfer coefficient for Plant 3 

Since the results of method 1 were very poor, only method 2 is studied for variation of 
mass transfer coefficient. Variations of 77 and capacity mass transfer coefficient {K') for 
Plant 3 can be observed from Fig. 5.42 (c) and Fig. 5.43 (d), respectively. It can be seen 
that, similar to plants 1 and 2 , K almost linearly decreases with time. 
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Graphs are drawn for heats with nearly same initial and 
final carbon vlaues 

(a) Co=0.369 Ct=0.044 (b) Co=0.369 Ct=0.044 

(c) Co=0.365 Ct=0.043 (d) Co=0.369 Ct=0.046 

Figure 5.42: Variation of r? with time for method 2 of GA-decarb models. 
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Capacity mass transfer coefficient (K,‘per sec) 


0.024 


0.022 


0.02 


0.018 


0.016 


0.014 


0.012 


0.01 

0 20 40 60 80 100 120 140 160 180 

Time (s) 

Graphs are drawn for heats with nearly same initial and 
final carbon values 

(a) Co = 0.369 and Ct=0.044 (b) Co = 0.369 and Ct=0.044 
(c) Co = 0.365 and Ct=0.043 (d) Co = 0.369 and Ct=0.046 

Figure 5.43: Variation of capacity mass transfer coefficient, K' with time for method 2 of 
GA-decarb models. 


(a) Plant 1, Case study 1(273 heats) 

(b) Plant 1, Case study 2 (156 heats) 

(c) Plant 2 (340 heats) - 

(d) Plant 3 (300 heats) 

... (d) 


- -...(c) 


(b) 


(a) 
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5.6.4 Comparision of capacity mass transfer coefficients for Plant 1, 
Plant 2 and Plant 3 

It can be observed from model 5.4 (method 1 , Case Study 2) of Plant 1 and model 5.4 
(method 1) of Plant 2 that cni and rj values ai'e larger in the latter model with 0.0084 and 
0.291, respectively. Because of the larger values of ai and rj, K' has greater magnitude 
as shown in Fig. 5.41 (c). This indicates that agitation of bath is higher in case of Plant 2 
(200 tn converter) than in Plant 1 (300 tn converter). This may be due to low lance height 
or less number of nozzles. 

After about 100 s, 77 reaches to almost same value for Plant 1, Plant 2 and Plant 3 
as shown in Fig. 5.42. The mass transfer coefficient in Fig. 5.43 (a) and (b) decreases 
very slowly with time. Hence, it can be concluded that lance height variation is properly 
adjusted in these cases. But in case of Fig. 5.43 (c) and (d), there is small amount of 
variation in K', indicating that lance height variation is not adjusted properly. This is 
specially true in the case of Plant 2 as in this case lance height is kept constant through out 
the blow. The values of K' and rj at time t=40s are shown in Table 5.20.The higher value 
of K'for Plant 3 is because of larger values of ai and 0:2 (Table 5.19) with 0.0195 and 
0.0175, respectively. Hence, it is observed from Fig. 5.43 that K' increases from Plant 
1 to Plant 3 (decreasing converter capacity). This indicates that in a 100 ton converter of 
Plant 3 bath is more agitated than in 200 ton converter of Plant 2 and 300 ton converter of 
Plant 1. This aspect is clear from Fig. 5.42. From this ftigure it can also be observed that 
there is a sharp drop in rj value in of Plant 2 and Plant 3 where as in Plant 1 the drop is 
more steady. This may be partially attributed to the changes in foaming nature of slag in 
the initial stages. After about 120 s, r? attains nearly the same value in Plant 1, Plant 2 and 
Plant 3 (Fig. 5.42). From Fig. 5.43, the value of K' is highest in the smallest converter 
(100 1 ) and smallest in 300 1 converter. 
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Table 5.19: Comparison of parameters in capacity mass trasfer coefficient 


"■Method 

Parameter 

Plant 1 Plant 2 

Plant 3 



Case study 1 

Case study 2 


1 

V 

0.2267 

0.2477 0.2524 

- 


ai 

0.00677 

0.0065 0.009 

- 

2 

ai 

0.00677 

0.0065 0.017 

0.0195 


Oi2 

0.00063 

0.00075 0.0014 

0.00175 


Table 5.20: Capacity mass transfer coefficient (/ir')and 77 at time t=40s. 


Parameter 

Plant 1 

Plant 2 

Plant 3 


Case Study 1 

Case Study 2 



V 

0.0092 

0.011 

0.015 

0.011 

K' 

0.01136 

0.01138 

0.0176 

0.0198 
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Chapter 6 

Conclusions and Suggestions for 
Further work 

6.1 Conclusions 

1. Models are developed for prediction of end point carbon, temperature, phosphorus, 
manganese and dissolved oxygen for three plants Plant 1, Plant 2 and Plant 3. Best models 
for end point carbon prediction are as summarized below and corresponding results are 
shown in appendix 2, appendix 3 and appendix 4 


Plant Model 

Model 

Parameters 

Statistics 

Remarks 

No 

equation 

used 



1 5.3 

d^ = iMo + 

Orel 

m=0.5796 

Case 

(reg) 


022 

c=0.0245 

study 2 



hlans2 

R=0.76 (0.7617) 

appendix : 



dolol 

(7=0.004 (0.0041) 




Co 

F=212.882 
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Plant 


2 


3 


Model 

Model 

Parameters 

Statistics 

No 

equation 

used 


5.3 

C't = 1^0 + MlCt+ 

0222, c2 

m=0.5498 

(reg) 

Er="2 

022 * Co 

c=0.0185 



reflux 

R=0.74 (0.7414) 



022, Co 

a-=0.005 (0.0055) 




F=389.423 

5.4 

dt = Ho + HiCt 

basicity 

m= 0.1929 

(reg, 

+ 

To, 022 

c= 0.0410 

met- 



R=0.44 (0.4374) 

hod 2) 



<r=0.004 (0.0038) 




F=70.494 


Remarks 


appendix 3 


appendix 4 


Best models for temperature prediction are as summarized below 


Plant Model 

Model 

Parameters 

Statistics 

Remarks 

No 

equation 

used 



1 5.6 

T = /Xo+ 

hot ratio 

m=0.8146 

Case 



022, To 

c=307.76 

study 2 



dolo2 

R=0.90 (0.9027) 

appendix 2 



slag2 

It 

b 





F=677.6145 
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Plant Model 

Model 

Parameters 

Statistics 

Remarks 

No 

equation 

used 



2 5.6 

T = /io+ 

022 

m=0.5269 

appendix 3 


EEiV-X 

To 

c=786.211 




Co 

R=0.72 (0.7241) 




flux 

II 

OO 





F=372.546 


3 5.6 

T = /ioH- 

022 

m=0.5351 

appendix 4 



To 

c=785.333 




raw dolo2 

R=0.73 (0.7310) 




hlans2 

a=9 





F=342.074 



Best models for phosphorus prediction are as summarized below 


Plant 

Model 

Model 

Parameters 

Statistics 

Remarks 


No 

equation 

used 



1 

5.11 

Pt ~ P>lPo 

SVOl.Toct 

m=0.8792 

Case 



+ Yh= 2 Pi ■ 

htr *p 

c=0.0014 

study 2 




slg2,Mn2 

R=0.93 (0.9253) 

appendix 2 




hlans2 

(7=0.0007 





Po 

F=917.438 


3 

5.11 

= Po + PlPo 

svol,llife 

m=0.6823 

appendix 3 




Co, To 

c=0.0024 





022, Mn2 

R=0.81 (0.8113) 





Tact 

(7=0.001 





Po 

F=574.0 
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Best models for manganese prediction are as sununarized below 


Plant Model 

Model 

Parameters 

Statistics 

Remarks 

No 

equation 

used 



1 5.14 

Mnj = /io + Z^iMno 

Pt,Co 

m=0.9237 

Case 


+ Z)r=2^ l^i • 

022,Mno 

c=0.0115 

study 2 



hlans2 

R=0.96 (0.9603) 

appendix 2 



dolo2 

(7=0.0065 





F=1827.27 


3 5.13 

Mn[ = //o + jJ'iMno 

raw dolo2 

m=0.6848 

appendix 3 


+ Y1Z=2 Ml ■ 

Mrio 

c=0.0456 




022 

R=0.83 (0.8274) 




hlans2 

(7=0.012 





F=647 



Best models for dissolved oxygen prediction are as summarized below 


Plant Model 

Model 

Parameters 

Statistics 

Remarks 

No 

equation 

used 



1 5.8 

[0] = /Lio+ 

Basicity 

m=0.761 

Case 


Y:=lP^^x^ 

022, Co 

c=101.255 

study 2 



To, dolo2 

R=0.87 (0.8724) 

appendix 2 



hlans2 

(7=34 





F=490.34061 
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Plant Model 

Model 

Parameters 

Statistics 

Remarks 

No 

equation 

used 



2 5.8 

[0] = fJ,o + 

Tact 

m=0.6695 

appendix 3 



To 

c=296.0932 




Cq 

R=0.82 (0.8182) 




022 

a=140 





F=684.434 


3 5.8 

[0] = /Xo+ 

basicity 

m=0.3571 

appendix 4 



To 

c=387.1967 




raw dolo2 

R=0.60 (0.5978) 




llife 

a=n 




Tact 

F=165.1967 



2. If the data for phosphorus and manganese are available at sublance point then the best 
model for prediction of end point carbon is model 5.3 (a) (reg) given by 


c't = - 8.219 * 10 “^ * Orel + 1.869 * 10 "“^ * + 114.642 * 10 ”® * 

84.317 * 10 "^ * Muo - 1.01 *po- 2.498 * 10 "^ * 022 - 0.2346 

with R=0.82, a =0.004%, F=308.958. This emphasises the point that during second blow 
carbon, manganese and phosphoms are simultaniosly removed and hence the variables 
phosphoms and manganese affect oxygen distribution significantly which, in turn, affects 
the kinetics of oxidation of carbon as well. 

3. It is concluded from GA-decarb models that capacity mass transfer coefficient depends 
on stirring energy and the value of rj decreases as a function of time in the range 0.005 to 
0.032. However, by properly adjusting lance height, mass transfer can be controlled such 
that the variation is reduced and becomes nearly constant. 

4. Capacity mass transfer coefficient ranges from 0.01 16408- 0.01 10834 for Plant 1 , 
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0.0188207 - 0.0170647 for Plant 2 and 0.0208267-0.0195667 (l/s)for Plant 3. The 
capacity mass transfer coefficient is highest for the smallest converter Plant 3 and least 
for largest converter Plant 1 . 

6. It can be observed from GA-decarb models that capacity mass transfer coefficient 
decreases with time very slowly as shown in method 2, indicating fiiat bath agitation is 
somewhat steady. More over GA-decarb models did not take into account the oxygen 
consumed for phosphorus and manganese removal .-Hence, linear models are performing 
better. 
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6.2 Suggestions for further work 


] . The effect of lance height vanaiion can be studied by incorporating the variation of 
effective gas flow rate (Qgf^-) and this effect can be included in GA-decarb models. 

2. A dynamic model can be made il waste gas analysis is available as a function of time. 

3. These models can be further tuned by using artificial intelligence techniques such as 
neural nets and fuzzy logic. 

4. The effect of lance nozzle design and blowing patterns for same converter need to be 
studied. 

5. Effect of droplet decarburization on slag foaming and post combustion can be 
incorporated. 

6. The dynamic effects of coolants added during the blow can be incorporated. 

7. A window based control model can be developed. 

8. Effects of heat loss and lining life and wear pattern can be incorporated by time series 
analysis, if data for sequential heats is available. 
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Appendix 1 

Specific stirring energy E% (WW ^ can be calculated from the expression given below. 


E'' 


V 


In 




i+Ml 




( 1 ) 


Where, 


Qb = gas flow rate = 3 m^/min, 

Ti = temperature of liquid metal 
p = density of metal = 7000 kg /m^, 

V = volume of liquid metal = (300 . 1000) / 7000 = 42.86 m^ 
G = acceleration due to gravity = 9.81 m/s^ 

H = bath depth = 1.8 m, 

P = ambient pressure = 1.0135 x 10^ Pa, 

Tm = temperature of gas = 298 °K. 


Upon inserting the values of different parameters in Eq.(l) we get. 


r 6. 18x3x1 873 

I. 

("7000x9.81x1.80 \ 

1 

ft 298 Y 


ln<l + 


> + 

1 

^ 42.86 , 

1 

[ 1.0135x10^ ] 

1 

1 1823 J_ 


= 884.87W /m^ 

Discussion:- It is assumed in the above calculation that bath depth does not 
change during bottom stirring owing to finite residence time of gas bubbles in 
liquid metal .Actually bath depth (or height) increases or bath level rises due to 
gas bubbles in metal ( see Example 4).Further , according to Eq.(l) the specific 
stirring energy is directly proportional to gas flow rate Qb.It also increases with 
bath height though not in direct proportion ,for example if the bath height is 
decreased to 1.42 m then corresponding stirring energy decreases to 818.27 W/m 
only. 


The total flow rate of oxygen is given by 
= 976.53 mV min. 



JMow the expression for specific energy received by metal from top of jet is , 

E“ = [ 6.32 X 10 cos O . M ] / V n^ dt LH (2) 

Where. 

0 = 14° , Q, = 976.53 mVmin , M = 32, V = 42.86 mVmin, n = 5 ,dt = 0.0466 m, 
lance height LH is 2 m, 

Substituting these values in to Eq(2), 

E°, = 84265.16 Watts /m^ 

Since only 10 % of the energy of jet is absorbed by bath the effective energy 
received is 8426.52 W/m^. 

Discussion:- Specific stirring energy received from top jet .varies inversely with 
lance height and square of number of nozzles; thus four hole lance will deliver 
approximately 1.56 times more energy than a five hole lance. Due to this reason , 
to properly utilize the jet energy, lance height is decreased as number of nozzles 
increase. Further as the number of nozzles increases the chances of jet 
interference increases. To avoid this nozzle angle (0) is increased as the number 
of nozzles increases. Thus, nozzle design and blowing regime are inter- 
dependent. 


_.(Kc)-=- 8.9X10-^ (Qeff)°^ 
where total stirring energy is 
Qeff = 1227.4 + 8552.54 watts/m^ 

(Kc) = 8.9 XIO"^ x( 9779.94)°^ 

= 0.014 

Where Kc = 0.014 m^/s (volumetric mass transfer coefficient). 
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Actual and Predicted end point parameters of Plant 3 
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